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SECTION  I 


INTRODUCTION 


At  present  there  is  an  intense  interest  among  many  re- 
search groups  in  learning  more  about  charges  and  currents 
induced  on  metallic  surfaces  of  objects  by  incident  electro- 
magnetic waves.  One  object  of  major  interest  is  the  aircraft. 
Of  course,  even  the  simplest  aircraft  is  an  extremely  complex 
body  when  viewed  in  the  context  of  a scatterer  in  an  electro- 
magnetic boundary  value  problem.  In  an  effort  to  obtain  at 
least  a cursory  understanding  of  the  behavior  of  charge  and 
current  induced  on  an  aircraft  surface,  several  research 
workers  have  replaced  the  actual  aircraft  by  a wire  model 
configured  after  the  general  shape  of  an  aircraft.  Then  the 
current  on  the  model  is  calculated  by  means  of  presently 
available  techniques,  with  the  hope  that  the  current  and 
charge  on  the  model  bear  a global  semblance  of  those  on  the 
actual  aircraft. 

Numerous  methods  [1-8]  for  calculating  currents  and 
charges  induced  on  thin-wire  structures  by  incident  electro- 
magnetic waves  have  been  developed  in  recent  years.  These 
methods  are  all,  of  course,  based  upon  Maxwell's  equations 
but  they  incorporate  a variety  of  approximations  to  achieve 
results.  Obviously,  the  exact  boundary  conditions  of 
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electromagnetic  theory  are  well  known  but,  in  dealing  with 
complex  thin-wire  structures,  one  employs  appropriate  simpli- 
fications to  reduce  his  problem  to  one  which  is  practically 
tractable.  For  a wire  ensemble  with  no  intersecting  wires 
(no  junctions),  the  traditional  assumptions  of  thin-wire 
theory  suffice  and  lead  to  accurate  solutions.  But  the 
question  of  what  boundary  conditions  to  impose  to  account  for 
the  confluence  of  two  or  more  wires  remains  partially  un- 
answered. The  conditions  which  are  imposed  by  present 
researchers  are  not  mutually  consistent. 

Since  all  junction  conditions  suggested  to  date  are 
more  or  less  approximate,  the  goals  that  one  seeks  to  achieve 
by  the  application  of  a given  junction  condition  should  be 
delineated.  Aside  from  ease  of  implementation  and  incor- 
poration into  the  solution  method,  the  conditions  should  at 
least  ensure  accurate  axial  linear  charge  and  current  densities 
on  all  parts  of  a given  structure  which  are  remote,  in  terms 
of  numbers  of  wire  radii,  from  junctions  of  wires.  With  this 
met,  possibly  the  next  most  important  objective  is  to  satisfy 
the  electromagnetic  boundary  conditions  in  some  sense  in  the 
junction  region.  At  the  present  time  it  is  not  known  to  what 
extent  these  two  goals  can  be  achieved  simultaneously  nor  is 
it  known  whether  or  not  exact  conditions  can  be  found. 

This  report  describes  a research  project  which  is  an 
initial  investigation  of  the  junction  conditions.  Within 
thin-wire  limitation,  the  correct  approximate  boundary 


conditions  are  developed  and  are  subjected  to  various  theo- 
retical and  experimental  test,  but  this  work  should  not  be 
looked  upon  as  final.  More  experimental  verification  is 
highly  desirous. 

The  purpose  of  the  study  described  in  this  report  has 
been  to  investigate  the  thin-wire  junction  conditions  and 
the  distributions  of  current  and  charge  on  wire  structures. 
Included  in  the  report  are  analytical,  numerical,  and 
experimental  investigations  of  the  behavior  of  current  and 
charge  on  a wire  with  a discontinuous  radius--the  so-called 
s tepped-radius  problem.  Measured  currents  and  charges  on 
wire  structures  with  bends  and  junctions  are  given,  and  an 
iterative  analysis  of  the  bent-wire  scatterer  is  included. 

In  Section  II  is  described  an  analytical  investigation 
of  the  fundamental  junction  created  by  the  confluence  of 
two  coaxial  wires  of  different  radii.  Based  on  theoretical 
considerations  of  this  model,  the  two  correct  conditions 
which  must  obtain  at  such  a junction  are  developed.  In 
Section  III,  a numerical  procedure  is  presented  for  analysis 
of  the  stepped-radius  wire  both  as  a scatterer  and  as  an 
antenna.  Experimental  results  of  Section  VI  serve  to  support, 
in  a preliminary  fashion,  both  the  junction  conditions  set 
forth  in  Section  II  and  the  numerical  results  of  Section  III. 

Section  V is  devoted  to  a complete  description  of  the 
measurement  facility  and  apparatus  as  well  as  to  pertinent 
construction  details  of  experimental  models.  As  a vehicle 
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in  which  to  apply  the  junction  conditions  proffered  in  II, 
an  iterative  solution  procedure  is  developed  in  Section  IV 
for  the  bent-wire  scatterer.  With  the  above-mentioned 
conditions  incorporated  in  the  analysis,  the  calculated  zero- 
order  current  distributions  on  the  bent  wire  are  shown  to 
be  almost  identical  in  shape  with  measured  currents.  Not 
only  does  the  iterative  solution  provide  a means  for  ..ending 
additional  credence  to  the  boundary  conditions  but,  also, 
of  importance  as  well,  it  is  an  analytical  technique  from 
which  one  can  determine  the  correct  current  distributions  on 
the  bent-wire  in  a direct  and  simple  manner. 

Section  VII  is  a thorough  report  of  first  measurements 
of  current  and  charge  on  a bent  wire  illuminated  by  a plane 
electromagnetic  field.  Also  in  this  section  are  found  first 
measurements  on  a rather  complex  structure  having  a junction 
formed  by  the  joining  of  four  non-perpendicular  wires.  This 
latter  structure,  dubbed  the  V-cross  scatterer  by  the 
investigators,  is  intended  to  be  a crude  but  useful  model 
of  a swept-wing  aircraft. 

Lastly,  to  assess  the  influence  of  junction  geometry 
on  currents  and  charges  remote  from  the  junction,  measure- 
ments are  given  for  the  case  of  a bent  wire  with  a small 
sphere  at  the  bend.  These  data  show  the  effect  which  a 
protuberance,  e.g.,  engine  pod,  tank,  on  an  object  might 
cause.  Also,  they  are  somewhat  indicative  of  what  factors 
one  must  consider  in  a complete  electromagnetic  representation 
of  a prescribed  junction  geometry. 
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SECTION  II 


JUNCTION  CONDITIONS  FOR  ELECTRICALLY  THIN  CONDUCTORS 

1.  INTRODUCTION 

An  important  step  in  the  determination  of  the  distribution 
of  current  at  all  points  on  the  surfaces  of  metal  radiating 
structures  that  include  confluent  or  intersecting  wires  is 
the  correct  specification  of  the  currents  and  their  derivatives 
at  the  junctions.  Examples  of  such  junctions  are  (a)  the  point 
of  confluence  of  two  straight  conductors  at  an  arbitrary  angle 
including  the  straight  angle  when  the  conductors  have  different 
radii;  (b)  the  junction  of  a vertical  antenna  with  the  radial 
wires  of  a top  load;  (c)  the  point  of  intersection  of  two  crossed 
antennas  that  make  an  arbitrary  angle;  (d)  the  junctions  in  a 
f ish-b one- 1 ike  configuration  of  conductors;  (e)  the  apex  of  an 
umbrella  antenna.  The  antennas  that  include  junctions  may  be 
driven  elements  or  parasitic  antennas  that  act  as  receivers  or 
scatterers.  Of  particular  interest  are  the  parasitic  crossed 
conductors  that  have  been  used  to  simulate  an  aircraft  exposed 
to  a periodic  or  transient  electromagnetic  field  of  high  inten- 
sity. An  accurate  knowledge  of  the  surface  densities  of  current 
and  charge  on  the  conductors  permits  the  specification  of  the 
tangential  magnetic  and  the  normal  electric  fields.  These  are 
needed  in  the  evaluation  of  the  fields  that  penetrate  into  the 
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interior  of  an  aircraft  through  an  imperfectly  conducting  metal 
skin  or  through  cracks  and  apertures. 

The  conditions  imposed  in  the  literature  on  the  currents 
at  the  junctions  of  electrically  thin  conductors  are  numerous 
and  varied.  There  is  a general  agreement  that  Kirchhoff's 
current  law  on  the  sum  of  the  inwardly  directed  currents  must 
be  enforced,  that  is 


n 


k-1 


This  follows  from  the  fact  that  for  electrically  thin  conductors 
(with  ka^<<l  for  all  radii  a^  of  the  confluent  wires)  the  charge- 
able surfaces  belonging  specifically  to  the  junction--as  a 

region  distinct  from  the  wires--are  negligibly  small  since  they 

2 

are  of  the  order  of  (ka)  . As  a consequence,  the  conservation 
of  electric  charge  requires  the  total  current  into  the  junction 
to  be  zero.  But  (1)  provides  only  one  of  the  n conditions 
required  if  there  are  n conductors  that  meet  at  a junction. 

The  n-1  additional  conditions  that  are  needed  to  specify  the 
n currents  are  given  by  Miller  et  al . [1]  and  Tesche  [2]  to  be: 


where  the  prime  denotes  differentiation  in  the  direction  of  the 
wire  and  a^  is  the  radius  of  wire  i.  These  conditions  are 
justified  by  their  proponents  with  the  statement  that  the  "sur- 
face densities  of  charge  must  approach  the  same  value  as  an 
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observer  moves  along  a wire  towards  the  junction  regardless  of 
which  wire  he  is  on."  Note  that  the  equation  of  continuity  for 
conductor  i is, 


“ -jw2ira^r|^  (3) 

where  is  the  charge  per  unit  length,  the  charge  per  unit 
surface  area.  It  follows  that  (2)  is  equivalent  to: 


n-,  = 


= n. 


(4) 


where  r)^  is  the  surface  density  of  charge  on  conductor  i at  the 
junction.  Note  also  that  approximate  rotational  symmetry  about 
each  conductor  has  been  assumed.  Other  writers  [3-6]  impose 
no  condition  on  the  derivatives  of  the  currents  at  the  junction 
but  state  that  continuity  of  scalar  potential  along  the  surfaces, 
of  all  conductors  at  the  junction  is  to  be  used  instead  of  (1) . 
This  is  carried  out  explicitly  in  the  work  of  Butler  [5],  whereas, 
as  pointed  out  by  Logan  [7],  continuity  of  scalar  potential  is 
"inherent  in  the  integral  formulation  of  Chao  and  Strait"[6]. 

King  and  Wu  [8] note  that  the  scalar  potential  is  in  any  case 
continuous  everywhere  except  across  a double  layer  as  in  a 
delta-function  generator.  They  impose  the  conditions: 


4 


^2 


^3 


= I 


(5) 


which  differ  from  (2)  whenever  the  aj^s  are  not  all  the  same. 
Specifically,  (5)  is  equivalent  to: 


qi  - q2  - q3 


(6) 
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I and  not  (4).  In  (6)  q.  is  the  charge  per  unit  length  on  conduc- 

i ^ 

> tor  i at  the  junction.  Actual  computations  by  Tesche  [2]  and 

I King  and  Wu  [8]  have  been  limited  to  the  junction  of  conductors 

!' 

I with  equal  radii,  in  which  case  (3)  is  equivalent  to  (2).  The 

i conditions  (5)  have  been  verified  exper:'.mentally  by  Burton  and 

King  [9]  but  only  for  conductors  with  equal  radii. 

In  a completely  analytical  approach  such  as  that  of  King 
and  Wu  [8]  a solution  cannot  be  obtained  without  the  explicit 
' application  of  n conditions  at  the  junction  of  n conductors, 

viz.,  (1)  and  (5).  When  numerical  techniques  such  as  the  method 
of  moments  are  used,  a rigorous  treatment  of  a junction  is  not 
possible  and  various  artifices  are  employed.  Some  of  these 
enforce  Kirchhoff's  law  (1)  explicitly,  others  do  not.  In 
general,  the  current  entering  the  junction  from  each  conductor 
is  obtained  by  extrapolation  and  the  extrapolated  currents  then 
satisfy  (1)  more  or  less  accurately.  Most  of  the  numerical 
procedures  actually  impose  no  explicit  condition  on  the  deri- 
vatives of  the  currents  as  the  junction  is  approached  so  that 
it  is  not  surprising  that  condition  (5)  is  usually  violated, 
j That  this  is  the  case  is  evident  from  many  of  the  graphs  of 

current  distribution  given  by  Logan  [7],  in  which  the  slopes  of 
; both  the  real  and  imaginary  parts  of  the  currents  in  the  con- 

! fluent  conductors  differ  greatly  at  the  junction.  Since  dis- 

I continuities  in  the  slopes  of  the  components  of  current  can 

occur  only  across  a delta-function  generator,  and  since  such 
generators  are  not  assumed  to  be  present,  the  local  distributions 
of  charge  per  unit  length  and  the  associated  radial  electric 
fields  cannot  be  correct  near  the  junction. 
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The  explicit  application  of  (1)  and  n-1  additional  conditions 
at  the  junction  of  n electrically  thin  conductors  is  required  in 
order  to  obtain  correct  solutions  for  the  currents  and  the  asso- 
ciated charges.  Evidence  will  be  presented  in  the  following 
sections  to  show  that  conditions  (5)  or  (6)  but  not  (2)  or  (3) 
should  be  applied  at  the  confluence  of  conductors  with  unequal 
radii.  Conditions  of  continuity  imposed  on  the  scalar  potential 
or  components  of  the  vector  potential  are  obviously  redundant 
and  no  substitute  for  (5). 

2.  CHANGE  IN  THE  RADIUS  OF  A COAXIAL  LINE 

The  distributions  of  current  and  charge  per  unit  length  on 
the  inner  conductor  of  a coaxial  line  resemble  those  along  a 
dipole  or  monopole  antenna  in  the  sense  that  both  are  approxi- 
mately sinusoidal  in  their  leading  terms.  This  suggests  that 
the  local  properties  of  the  current  and  charge  per  unit  length  ! 

at  and  near  a discontinuity  in  radius  should  be  comparable  at 
least  so  long  as  the  cross-sectional  dimensions  are  electrically 
small.  Such  a discontinuity  is  shown  in  Fig.  1.  Specifically, 
line  1 consists  of  an  inner  conductor  with  radius  a^  and  shield 
with  inner  radius  b in  the  range  of  negative  values  of  the 
coordinate  z;  line  2 extends  from  z=0  to  positive  values  of  z 
with  an  inner  conductor  with  radius  a2>a^^  and  the  same  shield. 

Thus,  there  is  an  annular  step  in  the  inner  conductor  at  z=0 
from  radius  a^  to  radius  a2  while  the  shield  continues  smoothly. 
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The  properties  of  the  line  with  this  change  in  radius  of 


the  inner  conductor  have  been  determined  by  variational  methods 
[10].  Each  line  is  characterized  by  the  conventional  parameters 
and  a corrective  lumped  susceptance  B connected  across  the  line 
at  the  junction,  z=0.  The  parameters  of  the  line  are; 


'i  ” in(b/a^) 


Jl.  = 
1 


y^Jln  (b  /a^) 
2tt 


(7a) 


(7b) 


^ci  “ y “ (?q/2tt)  £n(b/a^)  (7c) 

ki  = ojyrp:  = = kg  (7d) 

where  i=l,2.  For  a line  with  a^=0.5  cm,  a2=1.5  cm  and  b=5.5  cm, 
Cj^=23.2xl0  farad/m,  C2  = 15.2xl0  farad/m,  Jl^  = 4.8xl0  ^ henry/m 

i!.2“2.6xl0  ^ henry/m,  Z^^=143.9  ohms,  and  Z^2“78.0  ohms. 

The  general  formula  for  the  shunt  susceptance  B is  given 
by  Mar cuvitz  [ 10 ] . When  (a2~a^; << (b-a^) , which  is  of  interest 
here,  it  has  the  form: 


2A^ (a2~a^) ^ 
Z^.X(b-a.  ) 


(b-ai)  ^ ^2(b-a^)' 

(b-a,)  a, -a- 


With  B«uC  and  the  use  of  (7a)  and  (7c),  the  shunt  capacitance 
C at  z=0  is: 


C 


*=1 


•^1^' 


IT  (b  -a^) 


(b-a^^) 

(b-a2) 


£n 


m 


(8b) 


where  c^  is  the  capacitance  per  unit  length  of  line  1.  In  (8a, b) 


A 


1 


(b-a2)^  £n(b/a^) 

[£n(b/a2)]^ 


(9) 


A2  is  defined  in  general  on  page  230  of  Marcuvitz  [10].  For 
present  purposes  the  range  of  interest  is  given  by: 

(b/a^)> (b/a2)>>l  (10) 

In  this  range. 


b Jin  (b  /a^  ) 

A,<<1;  A = ± (11) 

a2[^.n(b/a2)  ] 


With  these  values  and  the  use  of  (10),  (8b)  becomes: 


^1  ^n  (b  /ai) 

^2/  a2[Jln(b/a2)]^ 


S-n  (b  /a, ) 
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When  b/aj^  and  b/a^  are  sufficiently  large,  the  leading  term  is: 

<=  ■ '1  i;) 

In  general  this  is  a very  small  capacitance  that  can  be  ne- 
glected in  a zero-order  approximation.  For  example,  with 
a^*0.5  cm  and  a2=l-5  cm,  C/c^  as  given  by  (12b)  has  the  value, 

C/c^  = 2.1  X 10  ^ m (12c) 

where  c^^  is  the  capacitance  per  meter  of  line  1. 

If  a standing  wave  is  maintained  on  the  line  with  a maxi- 
mum voltage  at  z^^O, 


V^(z)  = V^(0)cos  kz  , z<0  (13a) 

V^iz)  = V2(0)cos  kz  , z>0  (13b) 

With  the  condition, 

V^(0)  = V2(0)  (lA) 

and  the  basic  transmission-line  equation, 

q^(z)  = c^V^(z),  i = 1,2  (15) 


the  following  relations  are  found  to  obtain: 
q„(0)  c,  )in(b/a  ) . 

qj^(O)  la(b/a^)  in  3.67  ^ 

where  the  numerical  values  apply  to  the  case  of  a^=0.5  cm, 
a2“1.5  cm,  b=5.5  cm.  Thus,  in  zero-order  approximation,  the 
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charge  per  unit  length  is  discontinuous  across  the  junction. 

This  condition  is  illustrated  schematically  in  Figs,  la  and  2a 
where  the  charges  and  the  associated  radial  electric  fields 
change  their  values  abruptly  at  z=0.  Actually,  of  course,  the 
change  from  qj^(z)  to  q^Cz)  is  not  discontinuous  but  occurs 
over  very  short  distances  on  each  side  of  the  junction  with 
magnitudes  of  the  order  of  b.  In  this  range  the  capacitance 
per  unit  length  is  not  constant,  i.e.  c^(z)?^c^  and  c^(,z)^C2‘ 

Also,  the  charges  per  unit  length  vary  rapidly  but  continuously 
from  the  value  qj^(z)Aq^(0)  at  z = -5b  to  the  value  q2(z)=q,(0) 
at  z=5b.  Since  kb<<l,  it  is  evident  that 

q2(5b)  C2  iln(b/aj^) 

q^(-5b)  “ c^  inib/s.^) 

The  continuous  distributions  q^(z)  and  q2(z)  and  the  associ- 
ated electric  fields  are  shown  schematically  in  Figs,  lb, 2b. 

Note  that  outside  a narrow  range  on  each  side  of  z=0 , the 
distributions  of  q^^Cz)  and  q2(z)  are  the  same  as  in  Fig.  la. 
Inside  this  range  they  differ. 

In  both  Figs,  la  and  lb  the  small  corrective  capacitance 
C has  been  neglected.  If  it  were  included  in  Fig.  lb,  a small 
additional  charge  would  be  located  on  the  flat  annular  surface 
at  z=0  with  a consequent  very  small  shift  in  the  entire  standing- 
wave  pattern  toward  z=0.  The  actual  distributions  of  current 
and  charge  on  the  inner  and  outer  conductors  of  the  coaxial  line 
are  shown  schematically  in  Fig.  2b.  Note  that  q(z)  varies 
continuously  along  the  inner  surface  of  the  shield,  changing 
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OUTER  CONDUCTOR 


Figure  2.  (a)  Zero-Order,  (b)  Approximate  Actual  Charge  per  Unit 

Length  on  Inner  and  Outer  Conductors  of  Coaxial  Line  Near 
Change  in  Radius. 
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from  q^(z)=»c^V^(z)  to  q 2 ( z ) = C2  V2  ( z ) in  a very  short  distance 
on  each  side  of  the  junction  plane.  On  the  inner  conductor 
q(z)  drops  to  zero  at  the  sharp  inside  corner,  builds  up  along 
the  annular  surface  to  reach  a sharp  and  high  peak  at  the  out- 
side corner,  and  then  tapers  off  rapidly  to  reach  the  value 
q 2 ( z ) = C2^2 ( z ) . On  the  inside  surface  of  the  shield  the  change 
from  q^(z)  to  q2(z)  is  also  rapid  but  smooth.  The  magnitude 
of  the  total  positive  charge  in  the  range  -5b-z^5b  on  the  inner 
conductor  equals  that  of  the  negative  charge  on  the  outer 
conductor . 

The  significant  consequence  of  this  investigation  is  the 
condi t i on : 


q2(0)  in  b - Jin  a^ 
qj^(O)  ~ Jin  b - Jin  a2 


1;  Jin  b >>  Jin  a2 


(18) 


As  the  outer  conductor  recedes  from  the  inner  one  with  increas- 
ing b,  the  charges  per  unit  length  on  each  side  of  the  junction 
become  equal.  Although  this  relation  has  been  obtained  with 
a TEM  mode,  the  distributions  of  current  and  charge  per  unit 
length  along  an  antenna  are  sufficiently  similar  to  those  along 
the  inner  conductor  of  a coaxial  line,  that  it  seems  reasonable 
to  conclude  that 

q^(0)  = q2(0)  (19) 

for  the  antenna  with  a discontinuous  radius.  This  is,  of 
course,  a zero-order  representation  for  electrically  thin 
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conductors . 


3. 


CHANGE  IN  THE  CROSS  SECTION  OF  A SPHEROIDAL  ANTENNA 


Consider  a line  source  of  length  2h  that  extends  from 
z»-h  to  z=h  along  which  a suitable  distribution  of  generators 
maintains  a sinusoidal  current  of  the  form: 


I(z,t)  = I cosBz  coswt  ( 

’ m 

when  the  electrical  half  lengths  are  Sh=mr/2  with  n odd  and 


I(z,t)  = I sinBz  cosoJt  (I 

’ m 

when  n is  even.  The  associated  dis  t ributions  of  charges  f/er 
unit  length  are  obtained  with  the  equation  of  continuity 


9l (z , t)  9q (z  ,t)  _ „ 

3z  ' 3t 


They  are: 


q(z,t)  = (I  /c)  sinBz  sinwt  , n odd 
m 


q(z,t)  = (I  /c)  cosSz  sinojt,  n even  (221 

IQ 

g 

where  c=w/B=3xlO  m/sec  is  the  velocity  of  electromagnetic 
waves  in  air. 

The  exact  electromagnetic  field  at  all  points  generated 
by  these  distributions  of  current  and  charge  has  the  following 
components  in  the  cylindrical  coordinates  r , 1) , z for  n odd 
with  the  upper  sign,  n even  with  the  lower  signs 
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u I 

B^(t)  = - ± sin(cot  - BR2j^)] 


(23a) 


Ej.(t) 


- [^^^sin(o)t  - BR,,  ) ± ^^---^-  sin(03t  - BR-,  )]  (23b) 


4irr  R 


Ih 


Ih 


2h 


2h^ 


ih  2h 


(23c) 


.-7 


In  these  formulas  Uq=4ttx10  henry/m,  ^ ^ =yiI^Te^  = 1 2 0 it  ohms. 


Ih 


= 7(z-h)^  + , R =7(z+h)^  + 


(24) 


This  field  can  also  be  expressed  in  the  spheroidal  coor- 
dinates k ,k,  and  4>  with  the  ends  of  the  line  source  at  z = ±h 
e h 

as  the  foci  of  a family  of  prolate  spheroids  with  semi-major 
axes  a^=hk^  and  a family  of  orthogonal  hyperboloids  of  two 
sheets  with  s emi-con  j ugat  e axes  a^=hkj^.  These  are  defined  as 
follows  in  terms  of  the  distances  R^^  and  from  any  point 

P to  the  foci: 


^2h  ®lh  “ ^^e 


(25a) 


^2h  ■ ^Ih  ° ^^h 


(25b) 


The  spheroidal  coordinates 


ke  = ^e/^  = ; '^h  " 


(26) 


J 


are  illustrated  in  Fig.  3.  The  eccentricities  are  e^  for  the 
spheroids,  e^  for  the  hyperboloids.  The  coordinates  and 
have  the  following  ranges: 

1 < kg  ^ « ; -1  ^ kj^  ^ 1 (27) 

The  cylindrical  radial  coordinate  r can  be  expressed  in  terms 
of  the  spheroidal  coordinates  k^  and  kj^  as  follows: 

r = hyCk^  - 1) (1  - kj)  (28) 

With  (25a, b)  - (28),  the  cylindrical  components  in  (23a-c) 
can  be  transformed  into  the  following  spheroidal  components: 

cos(n7Tk  /2) 

' " 2^  ' r-r--  sin(u)t  - niTk^/2)  (29a) 

y(k^-l)d-kj) 

cos(nTTk,/2) 

E (k  ,k,  ,t)  = - , sin(wt  - nTTk  /2)  (29b) 

e e h (iTTn  „ Z ^ e 

7(k^-k^)(i-kj';) 

CqI  sin(mrk  /2) 

Ep(ke.kh.t)  = 2^  y=======  cos(o)t  - n7rk^/2)  (29c) 

y(k2-kj)(k2-l) 

for  n odd.  The  corresponding  expressions  for  n even  are: 
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(30a) 


Mq!^  sin(mTkj^/2) 
iTTh 


cos  (cot  - nTTk  /2 ) 


2irh 


cos  (cot 


nirk  /2 ) 
e 


(30b) 


C-I  cos  (niTk.  /2) 

-2^  ^ 


7(kf- 


sin(cot  - nTTk  /2)  (30c) 


In  these  formulas  E^(k^,kj^,t)  is  the  component  of  the  electric 

field  tangent  to  the  spheroid  defined  by  k “constant:  E (k  .k  t) 

e p e ’ h ’ 

is  the  component  of  the  electric  field  perpendicular  to  the 
spheroid  and,  hence,  tangent  to  the  hyperboloid  defined  by 

cons t an t . Graphs  showing  the  elliptically  polarized  electric 

field  for  n=l,2,  and  3 are  in  Figs.  8.4,  8.5, and  8.6  of  Theory 
of  Linear  Antennas  [11].  A detailed  description  of  the  field 
is  on  pages  540-546  of  [11].  Of  importance  here  are  the  facts 
that  B^(k^,kj^,t)  and  E^(k^,kj^,t)  are  always  and  everywhere  in 
phase  with  each  other  and  in  phase  quadrature  with  E ^ (k^  , k^^  , t ) . 
Each  spheroid,  k^^constant , is  a surface  of  constant  phase  or 
wave  front  that  expands  so  that  its  intersection  with  the  z 
axis  travels  with  the  velocity  of  light. 

Exactly  the  same  field  (29a-c)  or  (30a-c)  maintained 

by  a line  source,  k^=L,  with  the  current  (20a)  or  (20b)  can 

be  maintained  by  suitable  currents  on  a spheroid  k =k 

e el 
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for  all  values  of  greater  than  k^^.  The  required  currents 
and  the  associated  charges  can  be  obtained  by  requiring  them 
to  satisfy  the  boundary  conditions; 


fixB 


-POK 


(31a) 


ft  • E 


-n/e. 


(31b) 


where  K and  n are  the  surface  densities  of  current  and  charge, 
respectively,  on  the  surface  of  the  highly  conducting  spheroid 
defined  by  ^ externally  directed  normal  to  the 

Spheroidal  surface,  and  B and  E are  the  fields  with  components 
given  in  (29a-c)  or  (30a-c) . Specifically, 


Ig^^^el'^^h’^^  = 2iTr^K^(k^^,kj^,t)  = 2iTr  ^B  (k^^  ,k^^ , t ) /y^ 


(32) 


^e^^el ’^h ^ *^e  1 ’^h  ’ ^ ^ 0 


(33) 


whe  re 


(3A) 


is  the  radial  distance  from  the  axis  to  the  surface  of  the 
spheroid  at  each  value  of  k^ . When  (29a, c)  and  (30a, c) 

are  substituted  into  (32)  and  (33),  the  required  currents  on 
the  spheroidal  surface  are: 


n odd:  I (k  ,k,  ,t)  = -I  cos  (nirk.  /2 ) s in  (ut  - mrk  ./2)  (35a) 
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n even:  I (k  ,,k  ,t)  =■  -I  sin  (nTk,  /2)  cos  (ojt  - nirk  ,/2)  (35b) 

eeln  m n el 

The  associated  charges  per  unit  length  along  the  spheroidal 
surface  are: 


n odd : 


— ^ sin(mrk,  /2)  cos((jJt 

k^  -k^  ^ 

*^01  h 


mrkgj^/2)  (36a) 


n even: 


— ~ r cos(niTk  / 2)  sin(o)t 

el  h 


- n7rkg^/2)  (36b) 


A spheroidal  antenna  is  electrically  thin  when  its  semi- 
minor axis  b “hv/k^-l  satisfies  the  same  condition  as  the  radius 
of  a cylindrical  antenna.  That  is,  corresponding  to  the  con- 
dition 3a<<l  for  the  cylindrical  antenna  the  spheroidal  antenna 
must  obey  the  following  condition: 


6b  » 6hy  k^-1  = imr/2)  J k^-1  « 1 


(37a) 


This  is  equivalent  to; 


kf  - 1 5 0.1/n 
e 


k^  ^ 1 + 0.01/n  ; k^  ^ 1 + 0,005/n 


(37b) 
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It  is  evident  that  for  electrically  thin  spheroids  the  quantity 
2 

differs  very  little  from  1 so  that  the  factor  in  (36a, b) 


(38) 


is  also  near  1.  The  quantity  e(hj^)  is  shown  graphically  in 

Fig.  4.  It  is  seen  to  be  very  small  except  near  the  ends  of 

the  spheroid  where  kj^  = ±l.  The  trigonometric  functions 

s in (nTTkj^ /2 ) and  cos  (nTi'kj^/2)  are  not  able  to  take  full  account 

of  the  distributions  of  charge  as  the  ends  of  the  antenna  are 

approached.  The  quantity  >/(  1-kf ) / (k^-kf ) = l-e(k.  ) is  a 

n e (1  a 

weighting  factor  that  increases  with  decreasing  radius  of 
curvature  of  the  spheroid.  In  combination  with  the  trigonom- 
etric factors  it  provides  an  accurate  representation  of  the 
charge  density  on  the  spheroidal  antenna  including  the  high 
concentrations  at  the  ends. 

The  currents  and  charges  per  unit  length  on  a confocal 

spheroid  with  surface  defined  by  *^«»1  given  by  (35a, b) 

and  (36a, b)  with  k , replaced  by  k For  all  values  of  k >k  „ 

ei  e*  e e2 

the  fields  of  the  two  spheroids  with  different  k's  are  identical 

e 

The  amplitudes  of  the  currents  and  charges  on  the  two 

spheroidal  antennas  with  different  k's  are: 

e 


I 

e 


(kel«^h> 


-I 

m 


cos 

sin 


(mrkj^/2) 


(39a) 


Ie(k 


e2*^h^  “ 


-'m  s"^'^^kh/2) 


(39b) 


where  the  upper  trigonometric  function  is  for  n odd,  the  lower 
for  n even.  At  corresponding  values  of  on  the  two  spheroidal 

antennas  the  ratio  of  the  current  amplitudes  is 


^e  ^^e2  ’^h^ 


On  the  other  hand,  the  corresponding  ratio  of  the  charges  per 
unit  surface  length  is 


The  function  6(kj^)  is  also  shown  in  Fig.  4.  It  is  very  small 
except  near  the  ends  where  it  rises  rapidly  as  k^^  approaches  ±1. 
This  is  a consequence  of  the  fact  that  the  charge  per  unit  sur- 
face length  increases  most  quickly  where  the  radius  of  curvature 
of  the  thinner  spheroid  decreases  most  rapidly.  It  follows 
that  the  ratio  of  charges  per  unit  surface  length  on  the  thinner 
to  the  charges  on  the  thicker  antenna  rises  steeply  near  the 

ends.  At  the  center  of  the  antennas  (k,  *0)  the  ratio  of  the 

h 

charges  per  unit  surface  Idngth  is 
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1 . 00033 


(43) 


q 

q 


e 


e 


(k 


e 1 


(k 


e 2 


.0) 

ToT 


where  the  numerical  value  applies  to  spheroids  with  k^^^-l . 000166  , 

k „=1.0005.  At  the  ends  (k,=±l), 
e i h 


qe(ke2.il) 


1.732 


(44) 


which  is  significantly  greater  than  one.  The  rapidity  of  the 

rise  near  the  end  can  be  seen  from  the  fact  that  at  k, =0.995, 

n 


qe(kei. 0.995) 


0.990025 


0.990025 


1.065 


(45) 


Evidently  most  of  the  increase  occurs  within  0.005  of  the  ends. 

In  order  to  study  the  behavior  of  the  charge  near  a dis- 
continuity in  radius,  it  is  convenient  to  select  the  spheroid 
with  n=2,  i . e .,  Sh*TT,  for  which  the  maximum  of  charge  occurs 
at  the  center,  kj^  = 0.  The  currents  and  charges  are  given  by 
(39a, b)  and  (40a, b)  with  n = 2 for  two  spheroids  with  different 
k^s.  The  complete  electromagnetic  field  is  given  by  (30a-c) 
with  n=2;  it  is  shown  graphically  in  Fig.  8.5  of  [11].  The 

fields  of  the  two  antenn.s  are  identical  when  k >k  The 

e e2 

local  amplitude  of  E (k  ,k.)  near  k =0  is  shown  in  Fig.  5a 

p e h h 

where  the  upper  half  of  the  diagram  applies  to  the  thinner 

antenna  with  k =k  , , the  lower  half  to  the  thicker  antenna 
e el 
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Figure  5.  Charges  and  Electric  Fields  Near  Junction  of  Spheroidal 
Antennas  at  k,  = 0 where  E =0,  E Is  Maximum, 


with  k =k  Note  that  E (k  ,k,)=0  at  k,=0  and  is  very  small 

ee2  eeh  h 

near  For  clarity  the  values  of  k^^  and  k^^  actually 

used  in  the  diagram  are  much  greater  than  permitted  by  the 
condition  for  electrical  thinness.  The  electric  field  lines 
and  the  charges  from  which  they  emanate  are  shown  schemat- 
ically. The  charge  per  unit  surface  length,  “IgCkg.kj^),  is 
also  shown  along  each  surface. 

Although  Fig.  5a  was  drawn  to  illustrate  the  electric 
field  and  the  charges  near  the  point  of  maximum  for  one 
spheroid  with  surface  defined  by  in  the  upper  half  of  the 

diagram,  and  for  a second  spheroid  with  surface  defined  by 
k^2  io  the  lower  half,  the  representation  is  actually  a zero- 
order  approximation  of  the  fields  and  charges  of  a single 
spheroidal  antenna  with  surfaces  defined  by  different  values 
of  k^  in  the  upper  and  lower  half  spaces  and  a step  at  the 
junction  where  k^=0.  The  diagram  in  Fig.  5a  for  the  spher- 
oidal antenna  with  discontinuous  cross  section  at  its  center 
corresponds  to  Fig.  la  for  the  coaxial  line  with  a change 
in  radius  of  the  inner  conductor  at  z=0.  At  short  distances 
from  the  junction  the  electric  fields  are  unaffected  and  have 
the  form  and  magnitude  characteristics  of  an  unperturbed 
conductor.  Near  the  junction  the  charges  rearrange  them- 
selves near  the  step  in  the  sense  that  their  magnitudes  are 
reduced  at  the  inside  corner.  Increased  at  the  outside  corner 
with  appropriate  gradual  changes  along  the  adjacent  surfaces. 
These  are  Indicated  schematically  for  the  spheroid  in  Fig.  5b 
in  a manner  that  is  analogous  to  that  used  for  the  coaxial 
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Figure  6.  Charges  and  Electric  Fields  Near  Junction  of  Spheroidal 
Antennas  at  k,  = A where  E = 0,  E is  Maximum. 


line  in  Fig.  2b.  Note  that  at  short  distances  from  the  junc- 
tion the  fields  in  Figs.  5a  and  5b  are  the  same,  just  as  are 
those  in  Figs,  la  and  lb.  The  magnitude  of  the  charge  per 
unit  surface  length  is  shown  along  the  surfaces  in  Figs.  5a 
and  5b  in  a manner  similar  to  that  in  Figs.  2a  and  2b. 

The  important  observation  to  be  drawn  from  a study  of 
Figs.  5a  and  5b  is  that  the  ratio  of  the  zero-order  charges 
per  unit  surface  length  on  the  two  sides  of  the  junction  bet- 
ween electrically  thin  half  spheroids  with  different  eccen- 
tricities is  very  nearly  unity.  This  confirms  the  conclusions 
drawn  in  conjunction  with  (19)  for  a similar  junction  in  a 
coaxial  line. 


4.  REFINEMENT  OF  THE  THEORY 


The  discussion  of  the  junction  in  the  coaxial  line  and 
the  spheroid  is  highly  simplified  and  somewhat  unrealistic. 

In  the  coaxial  line  the  conclusions  depend  on  allowing  the 
radius  of  the  outer  conductor  to  approach  infinity  while 
preserving  the  purely  TEM  character  of  the  field.  This  implies 
that  the  cross  section  of  the  line  remains  electrically  small, 
so  that  the  limit  is  valid  only  at  zero  frequency.  Clearly, 
while  b should  become  large  compared  with  a^^  and  a^,  it  should 
not  exceed  a wavelength.  If  it  is  noted  that  (18)  is  equiv- 
alent to: 


!ci 

qi(0)  = 


(46) 
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and  if  it  is  recalled  that  in  coaxial  line  theory  the  char- 
acteristic impedance  plays  the  same  role  as  the  expansion 
parameter  4'=  2 [ 2,  n ( 2 / 6 a ) -0 . 5 7 7 2 ] in  antenna  theory,  the  indi- 
cation is  that  b should  approach  the  large  value  b=2/6  instead 
of  infinity.  In  this  case  (46)  becomes: 


q^CO)  £n(2/3a^) 

qj^(O)  ° ^ ° ZnCZ/Ba^) 


(47) 


If  the  radii  of  the  coaxial  lines  at  the  junction  are  made 
equal  to  the  semi-minor  axes  of  the  spheroids  used  in  (43), 
i . e .,  Sa2  = Sb^2~^  • ’ Baj^  = 6b^^  = 0 . 05  7 , (47)  gives: 


q,(0)  q^(0) 

qj^(O)  “ 


0.85 


(48) 


These  are  more  realistic  values  than  exactly  one. 

Throughout  the  discussion  of  the  spheroidal  antenna  it 
has  been  assumed  that  a suitable  excitation  exists  along  the 
antenna  to  maintain  currents  of  the  form  (39a, b)  with  the 
associated  charges  (40a, b).  Actually,  if  two  antennas  with 
different  semi-major  axes  and  ^g2“^^e2  same 

distances  between  foci,  viz.,  2h.  are  excited  by  identically 
distributed  fields  of  equal  amplitude  along  their  respective 
axes,  the  amplitudes  of  the  currents  will  not  be  the  same. 
Specifically,  as  shown  by  Ryder  [12], 

I (k  k ) 'v[Sh^k^.Y^(k^.)/Y'(k  )]  X (k  ) (49) 

n el  h el  n el  n el  n h 
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for  an  excitation  of  the  form 


n n 


1 - 


(50) 


In  these  formulas  X (k,  ) is  the  distribution  of  current  alone 

n h 

the  antenna  and  Y (k  , ) and  its  derivative  Y'(k  ,)  are  func- 

n e 1 n e 1 

tions  of  X (k,).  Y'(k  has  the  form 
n h n el 


k + 1 

Y ' (k  , ) = d Zn  T + c 

n e 1 k , - 1 

el 


(51) 


In  (51),  d and  c are  comparable  in  magnitude  in  general  whereas 
the  logarithm  is  quite  large  when  k^^  is  near  one.  However, 
d is  a function  of  frequency  and  goes  to  zero  at  the  reso- 
nances defined  by  Bh=mr/2  while  c remains  finite.  At  resonance, 

Y (k  ,)=exp(inTTk  /2),  X (k  ) =±  cos  (nTrk,  / 2 ) when  n is  odd, 
n e 1 e n n h 

X (k,  ) =±  s in  (nTTk,  / 2 ) when  n is  even.  It  follows  that  at  reso- 
n h h 

nance  and  when  exposed  to  identical  fields  of  the  form  in  (50) 
along  their  respective  axes,  the  amplitudes  I^  in  (39a, b) 
and  (40a, b)  for  the  antennas  must  be  replaced  by  I (k  ,) 
such  that  with  1^^  a constant. 


I(k.)=k-I^;I(k.,)=k_I,^ 
m el  elmO  m e2  e2m0 


(52) 


It  follows  further  that  (41)  and  (42)  must  be  replaced  by; 


‘el 


(53) 
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L 


In  particular,  the  spheroid  with  n = 2,  6h  = iT,  that  is  discon- 
tinuous at  k,=0  as  shown  in  Fig.  5,  will  have  zero-order 
h 

currents  given  by 


I (k  , ,k.  ) = -I^k  siniTk 
eelh  Oel  h 


(55a; 


I (k  o,k.  ) = -Ink  -sinTTk 
ee2h  0e2  n 


(55b; 


which  vanish  at  the  junction,  kj^  = 0.  The  corresponding 
charges  per  unit  length  are: 


1 - k: 


e 1 


el 


(56a 


1 - k 


k"„  - k 
e2 


COSTTk, 

2 h 


(56b 


where  c is  the  velocity  of  light.  At  the  junction,  k =0  and 


lg(k^.,0)  = ^e^^e2’^^  ~ ^ ’ '^e^^el’^^  ~ ~ ^0^^ 


(57) 


Thus,  in  the  resonant  antenna  with  the  junction  at  a charge 


maximum  and  current  zero,  with  each  half  identically  excited 


in  only  the  single  resonant  mode,  the  charges  per  unit  length 
at  the  junction  are  equal. 


With  the  same  excitation  but  not  at  resonance  the  loga- 
rithmic term  in  (51)  dominates  so  that  (52)  must  be  replaced 
by  : 

+ l)/(k^2  - 1)1  (58b) 

where  I is  a constant  and  Jin  [ k^+1 ) / (k^-1 ) ] corresponds  to 
the  expansion  parameter  '1'  of  a cylindrical  antenna.  It 
follows  that 


I^(k^^,kj^)  ~ 


’■e  ^^e2  ’’^h^ 


k^2  ^^el 


1)] 


(59) 


and 


Si  /“b2  - “h  - ^>1 


qe(‘^e2>‘^h> 


'e2 


- ‘^h  ^"[(‘^el 


(60) 


In  particular,  at  the  junction. 


9e^^e2’'^^  “^“f^^el 


(61) 


With  k^2  = l*0005  , k^  j^=  1 . 000 16  6 , (61)  gives: 
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r 


0.88 


(62) 


Note  that  this  value  is  in  good  agreement  with  (48)  for  an 
approximately  equivalent  cylindrical  conductor. 

When  the  spheroidal  antenna  is  not  excited  by  the  very 
specially  distributed  field  (50)  which  excites  only  the  n^^ 
mode,  but  is  immersed  in  an  incident  plane  wave,  the  total 
induced  current  is  the  sum  of  the  modes. 


I (k,  ) = Za  I (k  ) 
t h n n h 


(63) 


where  the  a^  are  amplitude  coefficients.  Even  when  the  fre- 
quency is  tuned  to  resonance  with  one  of  the  modes,  many 
higher  modes  will  also  be  present  but  with  much  smaller 
amplitudes.  It  is  to  be  expected,  therefore,  that  the  ideal 
conditions  leading  to  (53)  will  not  obtain  and  that  (53a, b) 
will  be  good  approximations.  This  means  that  the  condition 
(61)  on  the  charges  per  unit  length  at  a junction  should  be 
generally  useful. 

5.  CONCLUSION 


A careful  study  of  the  properties  of  currents  and  charges 
in  a coaxial  line  with  a step  in  its  inner  conductor  and  of  a 
spheroidal  antenna  with  a similar  step  leads  to  the  conclusion 
that  the  behavior  of  the  charges  per  unit  length  near  the 
junction  of  ttJO  cylindrical  conductors  with  different  radii 
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r rr- ^ 

is  not  simple  even  for  electrically  thin  conductors.  How- 
ever, a good  approximation  is  given  by  the  condition 

I 

'll^  = = = <^n^ 

where  4*^  is  the  expansion  parameter  for  conductor  i treated 
as  an  antenna.  The  definition  of  the  expansion  parameter  is 
somewhat  arbitrary,  but  is  always  a logarithmic  function. 

Useful  values  for  conductors  that  have  lengths  h that  satisfy 
6h-7r/2  are  given  by: 

T - 2[)ln(2/Ba)  - 0.5772  ] (65) 


For  shorter  antennas  the  appropriate  parameter  is 

4<  = 2 )ln(2h/a)  (66) 

where  a is  the  radius  of  the  conductor.  The  conditions  (64) 
and  the  Kirchhoff  condition  (1)  apply  to  junctions  of  elec- 
trically thin  conductors.  There  is  no  justification  for  a 
condition  of  the  form  (2)  or  (4).  The  conditions  (5)  or 
(6)  may  be  adequate  approximations  of  (64)  for  junctions  of 
conductors  that  are  sufficiently  thin  and  not  too  different 
in  cross-sectional  size. 

Since  in  thin-wire  theory  the  surface  area  of  a junction 
is  negligible,  and  each  conductor  is  treated  as  if  its  charges 
were  concentrated  at  an  average  position  on  the  axis,  the 
angle  at  which  two  conductors  meet  is  irrelevant.  This  is 
true  provided  the  diagonal  distance  across  a junction  is 
electrically  small,  and  not  only  the  diameter  of  the  thicker 
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SECTION  III 


ANALYSIS  OF  STEPPED-RADIUS  WIRE 
1.  INTRODUCTION 

In  this  section  Is  presented  an  analysis  of  the  stepped- 
radius  wire  structure  shown  in  Fig.  7.  This  structure,  with 
its  confluence  of  two  coaxial  wires  of  different  radii,  is 
perhaps  representative  of  the  simplest  thln-wire  configura- 
tion possessing  a junction.  For  this  reason,  the  behavior 
of  the  s tenped-r adius  wire  must  be  understood  before  more 
complicated  junction  problems  are  undertaken.  A numerical 
procedure  is  developed  for  solving  for  the  current  on  the 
s tepped-r ad ius  wire.  The  procedure  is  applied  to  both  a 
s tepped-radius  scatterer  and  an  antenna,  and  representative 
results  are  provided. 

In  the  present  investigation  of  the  stepped-radius  struc- 
ture of  Fig.  7 , the  wire  radii,  a and  b,  are  looked  upon  as 
being  very  small  relative  to  the  wavelength  X as  well  as  to 
the  wire  length  L.  These  restrictions,  common  in  thin-wire 
theorv,  assure  one  that  the  current  on  each  cylindrical  sur- 
face is  circumferentially  independent  and  that  it  may  be 
accounted  for  bv  the  total  axial  current, denoted  I on  the 

cvlinder  of  radius  a and  I,  on  that  of  radius  b.  These 

b 

quantities  are  shown  in  Fig.  7 , which  serves  to  define 
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geometric  quantities  of  interest. 

In  the  analysis  which  follows,  the  boundary  condition 
that  total  tangential  electric  field  be  zero  on  the  conducting 
surface  is  enforced  on  cylindrical  surfaces  but  not  on  the 
small  annular  surface  located  where  the  wire  radius  is  dis- 
continuous at  z=s.  This  annulus  is  a very  small  portion  of 
the  total  surface  area  of  the  s tepped-radius  wire  structure 
and,  therefore,  when  the  wires  satisfy  thin-wire  conditions, 
the  accuracy  of  a solution  for  current  should  not  suffer 
seriously  from  a failure  to  enforce  the  boundary  condition 
on  the  annulus.  Experience  with  thin-wire  analyses  supports 
this  simplification,  since  excellent  solution  accuracy  is 
obtainable  for  problems  involving  solid  wires  even  though 
the  boundary  condition  mentioned  above  is  seldom  enforced  on 
the  disks  at  the  wire  ends. 

A further  simplification  employed  here  is  that  one 
assumes  he  may  ignore  the  small  fraction  of  total  field  due 
to  the  induced  charge  and  current  which  reside  on  the  annulus. 
The  chargeable  annular  surface  is,  indeed,  small  compared  with 
the  cylinder  surfaces,  and  only  a relatively  small  fraction 
of  the  total  charge  on  the  structure  could  reside  on  the 
annulus.  Ignoring  the  charge  on  the  annulus  or,  equivalently, 
approximating  its  small  value  by  zero,  is  tantamount  to 
requiring  the  total  axial  current  to  be  continuous  at  the 
point  z=s  (junction)  where  the  two  wires  of  different  radii 
join.  If  the  axial  current  were  discontinuous  at  the  step, 
charge  would  be  deposited  on  the  annulus  in  order  that  the 


continuity  equation  be  honored.  Therefore,  under  the  assump- 
tions set  forth,  one  considers  only  the  induced  charge  and 
current  on  the  cylindrical  surfaces  and  enforces  the  boundary 

condition  E =0  only  on  these  two  surfaces, 
tan  ^ 

In  a scattering  problem,  one  requires  the  sum  of  the 

prescribed  incident  electric  field  and  the  scattered 
— 2 

electric  field  E , due  to  induced  sources,  to  have  zero 

tangential  component  evaluated  on  the  perfectly  conducting 

••  s ” i /N 

cylindrical  surfaces  of  the  structure:  (E  +E  ) *u^=0  on  the 

surfaces.  In  an  antenna  problem,  one  may  view  the  tangential 

incident  field  as  being  localized  at  the  emf  generator  and, 

again,  he  requires  this  specified  field  and  the  tangential 

scattered  field  to  add  to  zero  on  the  structure:  E^*u  +V6(z-z  )=0 

^ 8 

on  the  surfaces,  where  V6(z-z  ) is  a delta  function  of  strength 

8 

V at  due  to  a slice  generator  of  V volts  located  at  z^. 

Due  to  the  thin-wire  assumption  and  the  rotational  symmetry 
in  both  the  scattering  and  antenna  problems,  the  scattered 
electric  field  E is  independent  of  (j)  and  has  no  component  in 
the  0 direction.  Hence,  to  satisfy  the  boundary  condition 
that  the  total  tangential  field  be  zero  on  the  two  cylindrical 

g 

surfaces,  one  needs  only  to  calculate  E^(r,z),  the  axially 

S 1 

directed  scattered  field,  and  to  require  E^+E^=0  on  the  surfaces. 
From  basic  electromagnetic  theory  applied  to  the 

g 

structureof  Fig.  7,  one  may  readily  calculate  E^(r,z)  from 

E®  = -joiA  - (67) 

Z Z d Z 
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where  ( r , z ) is  the  z-component  of  the  magnetic  vector  poten- 
tial and  $ ( r , z ) is  the  electric  scalar  potential,  both  of  which 
are  readily  determined  from  the  induced  current  and  charge 
on  the  structure.  The  harmonic  time  variation  of  angular  fre- 
quency CO  is  suppressed  in  (67).  In  view  of  the  assumptions 
and  simplifications  discussed  above,  and  $ are  calculated 
from  the  current  and  charge,  respectively,  which  reside  wholly 
on  the  cylindrical  surfaces;  the  contributions  to  $ from 
charge  induced  on  other  surfaces  are  ignored.  The  boundary 
condition  discussed  above  now  can  be  written 


E® (a,z)  + E^(a,z) 


ze(s,L/2) 


(68a) 


and 


E®(b,z)  + E^(b,z)  = 0 , ze(-L/2,s) 
z z 

for  the  scattering  problem  and 


E®(a,z)  + V5(z-z  ) = 0 ze(s,L/2) 

z g 

and 

E®(b,z)  + V5(z-z  ) = 0 ze(-L/2,s) 

z g 


for  the  antenna  problem. 


(68b) 


(69a) 


(69b) 
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2.  NUMERICAL  SOLUTION  PROCEDURE 

From  Eqs.  (67)  and  (68)  or  (69),  with  and  $ written 
as  functions  of  current,  one  may  readily  obtain  an  integro- 
differential  equation  characterizing  the  current  on  the 
s t epD ed- r ad  i us  structure.  Then,  of  course,  this  integro- 
differential  equation  may  be  solved  for  the  unknown  current. 
The  solution  procedure  employed  here  is  essentially  the  moment 
method  [13]  but  with  a change  in  the  usual  sequence  of  steps 
of  this  method  which  enables  one  to  gain  interesting  insight 
into  the  nature  the  numerical  procedure.  The  alteration 
alluded  to  is  that  one  performs  the  testing  [13]  of  the  equa- 
tion to  be  solved  before  the  unknown  current  is  anproximated 
as  a linear  combination  of  the  basis  set.  The  testing  pro- 
cedure leads  to  a system  of  linear  equations  which  represents 
a discretized  equivalent  of  the  original  i n t e g r o-d  i f f e r en t i a 1 
equation.  The  numerical  solution  orocedure  is  described  below, 
in  general,  for  a representative  segment  of  wire  beginning  at 
and  extending  to  z^^,  and,  subsequently,  the  method  is 
aoolied  to  the  s t epp ed- r ad ius  structure, 
a . Test ing 

Pursuant  to  the  establishment  of  a linear  svstem  of 

eouations  from  which  one  raav  ultimately  obtain  a solution  to 

Eo  . (67),  one  equates  the  co  r r es  nond  i n g pro-jections  of  the 

two  sides  of  Eo.  (67)  onto  the  snace  spanned  bv  the  testins 

set.  In  other  words,  if  T is  an  element  of  the  testing  set. 

m 

one  equates  the  scalar  nroduct  of  T and  the  left-hand  side 

in 
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of  (67)  to  that  of  T and  the  right-hand  side: 

m 


/-juA  - T > = T > 

N-"  2 3z  m/  \z  m/ 


(70) 


m = 1 , 2 , , . . , M 


where  the  scalar  product  is  defined 

L/2 


<t,  s>  - j 


f g*d  z 
z=-L  / 2 


(71) 


Because  they  are  particularly  amenable  to  numerical  solutions 
of  (67)  and  lead  to  stable  linear  systems,  the  set  of  piece- 

wise  linear  functions  is  selected  for  testing;  the  m^^  element 

I 

A (z)  of  this  testing  set  is  defined 
m 


4 


0 

m 


- 


"^^^m-l'"m+l^ 


(72) 


where 


A 


( z. 


z 

m- 


1 


) 


(73) 


Implicit  in  the  above,  and  depicted  in  Fig.  8 . is  the  parti- 
tioning of  the  segment  of  wire  of  length  z^^-z^  into  M-1 
subintervals,  each  of  length  A=  ( Zj^- z ^ ) / (M- 1 ) . 
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Figure  8.  Testing  Triangle  Functions 


Since  the  testing  function  A'  is  zero  outside  the  sub- 

m 

domain  ze  ( z^_  ) , one  sees  immediately  that  (70),  subject 

j 

to  (71)  and  (72),  reduces  to  (T  =A') 

m m 


>1^-* 


^)(r,z^)  - ^ J*^(r  ,z)dz  -joj  ( r , z ) ( 2 ) d z 


2*z . 


2 = z . 


z 

I 


Z-Z  . 


E®(r,z)A^(z)dz  , (75a) 


(m=l) 


'm+1 


‘m+1 


J*<i>(r,z)dz  - j J’i'(r,z)dz  -jw  jA^(r,z)A^ 


( z ) d z 


z = z 


m- 1 


2 = Z 


Z = Z 


m- 1 


'm+i 


/ 


E^(r,z)A  (z)dz  , 
2 tn 


(75b) 


'm- 1 


m = 2,3,4,. 


and 


5 4 


r 


M 


-$(r 


,2m)  J^(r,z)dz  -ju)  J 


z = z 


M-1 


A^(r,z)AM(z)d2 


z = z 


M-1 


'M 


/ 


E(r,z)A(z)dz 
z M 


z = z 


M-1 


(75c) 


(m=M) 


The  testing  functions,  often  called  triangle  functions,  are 

illustrated  in  Fig.  8 . For  m= 2 , 3 , . . . , M- 1 , the  testing  functions 

are  interior  full  triangles  and  are  associated  with  Eq.  (75b). 

On  the  other  hand,  due  to  the  integration  limits  on  (71), 

£ 

testing  with  A^  is  equivalent  to  testing  with  a hal f - t r i ang 1 e 

£ 

on  the  lower  end  while  use  of  Am  yields  the  same  results  as 
would  a half-triangle  on  the  upper  end  (Fig.  8 ). 

The  vector  potential  due  to  currents  on  a wire  is 
reasonably  slowly  varying  with  z.  Therefore,  with  A suffi- 
ciently small  relative  to  wavelength,  the  approximations  below 
are  quite  good,  and  they  significantly  lessen  the  computa- 
tional complexity  of  the  present  analysis: 


/ 


A^(r,z)A^(z)dz 


j A^(r,z^)  , A<<X 


(76a) 
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(76b) 


'm+1 


/ 


A (r,z)A  (z)dz  = AA  (r,z  ) , A<<X 

z m z m 


'm- 1 


'M 


/ 


A (r , z)A  (z) dz 
z n 


2 


'M-1 


(76c) 


The  same  type  approximation  can  be  applied  to  the  right-hand 
sides  of  (75),  but  is  usually  unnecessary  because  the  indicated 
integration  can  almost  always  be  performed  analytically.  How- 
ever, in  the  interest  of  notational  convenience,  the  following 
definitions  are  employed: 


V=(r,z^) 


/ 


E®(r,z)A^(z)dz 


(77a) 


V (r,z  ) 
m 


m+1 


E“(r,z)A  (z)dz 
z m 


'm-1 


m=2 , 3 , . . . ,M-1 


(77b) 


1*1 

= J 


E®(r ,z)A^(z)dz 


'M-1 


(77c) 
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With  the  approximations  (76)  and  definitions  (77),  the 


three  equations  of  (75)  can  be  written  in  the  form 


i. 

'I’(r,Zj^)  ~ J' ^’(r,2)dz  -jw  j 


z = z . 


Az  (r  ,Zj^)  = V (r  ,z^) 


(m=l) 


(78a) 


m 


'm+1 


/$(r,z)dz  - j /<t(r,z)dz  -jojAA  (r,z  ) = V®(r,z  ) , (78b) 

4A  I z m m 


Z*2 


m- 1 


m 


m=2,3 M-1 


M 


i /«<'. 


■$(r,z,_j)  + j |$(r,z)dz  -joj  j A^(r,Zj^)  = V®(r,z^)  , (78c) 


z = z 


M-1 


(m=M) 


b.  Current  Approximation 

The  next  step  in  the  numerical  solution  procedure  is  to 
approximate  the  unknown  current  I(z)  in  a series  of  the  form 

N 

n=  1 
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I(z) 


(79) 


where  each  i (z)  is  a known  element  of  a basis  set  selected 
n 

for  representing  the  current,  and  I is  the  unknown  complex 

n 

t h 

coefficient  of  the  n term  in  the  series  (79).  It  is  the 

set  of  coefficients  {l  } that  one  seeks  to  determine  in  the 

n 

numerical  procedure.  and  are  written. .as  functions  of 

the  terms  in  (79)  and  then  are  substituted  into  (78),  which 
leads  to  a system  of  linear  equations  that,  subject  to  bound- 
ary conditions,  can  be  solved  for  {l  } . 

n 

A very  useful  basis  set  in  terms  of  which  to  represent 

the  current  on  the  s t epped- r ad ius  wire  is  that  comprising 

the  so-called  pulse  functions  (i  =p  ): 

n *^n 


P„(z) 


■ 


2£ (z  - A/2 , z + A/2) 
n n 


zi(z  - A/2 , z + A/2) 
n n 


(80) 


where  z^  locates  the  center  of  the  n pulse.  A represen- 
tative pulse  is  shown  in  Fig.  9 together  with  half-pulses 
on  the  ends  of  the  interval  (z^,z^^).  The  number  of  pulses, 
including  the  two  half-pulses,  is  chosen  to  be  equal  to  M, 
the  number  of  testing  triangles. 

c.  A and  $ in  Terms  of  I 
z n 

The  vector  potential  A^  due  to  the  approximation  given 
in  (79)  is  conveniently  written 

M 


A^(r  ,z) 


/ I A (r , z , R) 
n n 

n=  1 


(81) 
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'1 


In  (81)  A is  the  contribution  to  the  total  vector  potential 
n 

from  the  n*"^  pulse: 


z^+A/2 


Aj^(r,z,R)  = |g(z-z'  ,r,R)dz' 


/* 


(82a) 


z +A/2 
n 


A^(r , z ,R) 


/ 


g(z-z' ,r,R)dz'  , 


z ' = z - A / 2 
n 


n=2 , 3 , . , . ,M-1 


(82b) 


'M 


Aj^(r  , z ,R) 


4it 


/ 


g(z-z' ,r,R)dz' 


z '=z„-A/2 
M 


(82c) 


where  R is  the  radius  of  the  wire  on  which  resides  the  current 
I being  approximated  by  the  given  pulse  and  where  U is  the 
permeability  of  the  medium  in  which  the  wire  resides.  In 
(82)  the  kernel  is  given  by 


g(5 


■'•■x  - i / 


•jk[5^  + r-+R  — 2rR  cos(()'  ] 


({)  ' =-TT 


[C^+r^+R^-2rR  cos<p  ' ] ^ 


— <!$' 


(83) 


where  k is  2ir/>..  By  a suitable  change  of  variables  in  each 
expression  of  (82),  these  integrals  can  be  converted  to  the 
f o 1 lowing, wh i ch  are  far  better  suited  for  computation  than 
are  those  of  (82): 
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(r  , z ,R) 


g(z-z'-Zj^,r,R)dz' 


A (r , z , R) 
n 


— f 

4it  J 


Aj^(r,z,R) 


g(z-z'-z  ,r,R)dz' 
n 


z'*-A/2 


n = 2 , 3 , . . . ,M-1 


— f 

4tt  J 


g(z-z'-Zj_j,r,R)dz' 


z'=-A/2 


The  scalar  potential  <J>  due  to  the  charge  associated  with 
the  current  approximation  (79),  with  the  pulses  of  (80),  is 
to  be  calculated  as  a sum  of  partial  potentials. 


$(r,z)  = > I^$^(r,z,R)  , 


where,  of  course,  $ is  that  contribution  to  (J)  due  to  the 

n 

charge  related  to  p^.  As  suggested  in  Fig.  9 , current  on 
a cylinder  in  the  form  of  a pulse  "deposits"  a discrete  ring 
of  charge  at  the  upper  and  lower  ends  of  the  pulse,  where  it 
is  discontinuous.  Since  the  derivative  of  the  pulse  is  zero 
other  than  at  the  points  of  discontinuity,  no  other  charge, 
e.g.,  linear  charge  density,  is  associated  with  the  current 
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pulse.  The  total  charge  deposited  at  the  n pulse's  upper 


ead-point,  z^+t\/2,  is  - — and  that  at  its  lower  end-point, 

2.  -till,  is  + The  current  on  the  wire  surface  is  assumed 

N JO) 

to  be  uniform  around  the  wire  circumference  so  each  ring  of 


charge  is  of  constant  linear  density:  + j m 2 tt R 

Now,  each  partial  potential  <I>  is  seen  to  be 

n 


0^(r,z,R)  = 


j4^uie  [g(z-z^-A/2  ,r  ,R)  -g  (z-z^  , r , R)  j (86a) 


■f^(r,z,R)  = 


- ^g(z-z^-A/2,r,R)  -g  ( z - z^+A  / 2 , r , R ) , (86b) 


n=2,3, . . . ,M-1 


<l>j^(r,z,R) 


■ j4^  g(z-z^,r,R)  -g(z-z^+A/2,r,R)  (86c) 


where  e is  the  permittivity  of  the  medium  surrounding  the  wire 
structure . 

d.  Linear  System  of  Equations  for  S t epped-Rad ius  Scatterer 
If  of  (81)  and  <5  of  (85)  subject,  respectively,  to 
(82)  and  (86),  are  substituted  into  (78),  one  obtains  a system 
of  linear  algebraic  equations  having  the  current  coefficients 


I^'s  as  unknowns.  This  system  of  equations  comes  from  testing 
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(67)  with  the  triangle  testing  functions,  and  it  contains  the 
same  information  as  does  (67). 

At  this  point  we  adapt  the  system  of  equations  to  the 
s tepped-radius  wire  scatterer.  The  scheme  for  applying  the 
above  to  the  s tepped-radius  problem  is  suggested  graphically 
in  Fig.  10.  In  this  example,  the  portion  of  the  stepped- 
radius  wire  between  z=-L/2  and  z=s  (step  location)  is  divided 
into  six  subintervals  of  length  = [ s+L  / 2 ] / 6 , and  the  remaining 
portion  is  divided  into  three  subintervals  of  length 
A =(L/2-s)/3.  Next  (78)  is  applied  individually  to  the  in- 
terval  ze(-L/2,s)  and  to  ze(s,L/2),  and  the  currents  I.  and 

D 

are  approximated  by  pulses  in  the  re^oective  intervals. 

The  testing  triangles  and  the  pulses  are  observed  in  Fig. 10 
where  one  sees  that  the  coefficients  and  of  the  end 

half-pulses  are  made  equal  to  zero  so  that  Ij^(-L/2)=0  and 
I ^ ( L/ 2 ) =0-- cur r ent  boundary  conditions  at  free  wire  ends.  With 

3 

b 3 

Ij^  = I^Q  = 0,  the  number  of  unknowns  is  reduced  by  two,  and,  thus, 
(78a)  applied  in  (-L/2,s)  and  (78c)  applied  in  (s,L/2)  are 
deleted  to  reduce  correspondingly  the  number  of  equations  in 
the  linear  system.  Also,  notice  that  the  current  is  made 
continuous  at  the  step  by  the  composition  of  the  7^'^  pulse; 
this  pulse  occupies  an  interval  of  length  *i(A^+Aj^)  and  is  the 
combination  of  the  upper  half-pulse  on  the  1 ar ge r- r ad ius  wire 
and  the  lower  half-pulse  on  the  sma 1 1 e r- r ad i u s wire.  Com- 
bining the  two  half-pulses  at  the  step  into  i single  pulse 
(I-*I^»I)  reduces  the  number  of  unknowns  by  one,  and  a corre- 
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spending  reduction  is  achieved  in  the  linear  system  by  adding 
(78c),  applied  in  (-L/2,s),  to  (78a),  applied  in  (s,L/2). 


With  the  current  represented  by  pulses  and  with  the 
testing  of  (67)  (leading  to  (78))  suggested  by  the  triangles 
in  Fig.  10,  one  applies  (68)  to  in  (77)  and  enforces  (78) 
on  the  surface  of  the  larger-radius  wire  and  on  the  surface  of 
the  smal ler- r ad ius  wire.  The  linear  system  resulting  from 
the  above  applications  of  (78)  to  the  two  wires  in  the  struc- 
ture of  Fig.  10  can  be  readily  solved  for  the  current  coef- 
3 b 

ficients  I and  I . 

n n 

For  a general  s tepped-rad ius  scatterer,  one  would  divide 
(s+L/2)  into  segments  of  length  A^= (s+L/2 ) /N^  and  (L/2-s) 
into  segments  of  length  A^= ( L/ 2- s ) /N^ . Then,  applying 
the  scheme  described  above  for  the  example  of  Fig. 10, he  would 
establish  the  following  system  of  linear  algebraic  equations: 


(87a) 


(87b) 
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+ IZ^  + 
n mn  m 


n mn 


n = N, +2 

D 


"='Nb  + 2 N^+Nj^ 


where  and  are  the  weighting  coefficients  of  the  pulse 
approximations  for  I and  I,  , respectively.  The  elements  in 

cL  D 


the  above  linear  system  are  defined  in  (88). 


,z)A  (z)dz 
m 


^^(b.z)  [a^-s  + z] 


z = s-A, 


a 


(88b) 


r 


1 

i 

j 

X 
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For  a specified  incident  field  having  a z component  E^, 
the  system  of  equations  (87),  with  the  definitions  of  (88), 
can  be  solved  for  and  from  which  a pulse  approximation 

for  current  on  the  scatterer  is  readily  available, 
e.  Reduced  Kernel 

Although  the  exact  kernel  of  (83)  can  be  calculated  by 
means  of  elliptic  integrals  and  can  be  employed  in  the  com- 
putation of  quantities  above,  an  approximation  of  (83)  lessens 
the  computer  time  required  in  calculations. 

The  distance  from  a point  (r,0)  on  the  inner  circle  of 
Fig.  11  to  (R,c()')  on  the  outer  circle  is,  from  the  law  of 

U 

2 9 ^ 

cosines,  [r  +R~-2rR  cossj)']  . The  maximum  value  of  this  dis- 

2 2^  2 2 
tance  is  [r  +R  +2rR]  and  Che  minimum  is  [r  +R  -2rR]  from 

which  one  determines  the  median  distance  to  be 

Hence,  one  replaces  (r^+R^-2rR  cos(p')  in  (83)  by  (r‘"+R^)/2 
and  arrives  at  a reduced  or  average  kernel  for  the  stepped- 
radius  wire: 


2 21 
r +R^ 

2 . 


^-jk[C^+4(r2+R2) ] ' 

K(C,r,R)  = ^ ^ (89) 

[C  +*i(r"+R^)  ] 

The  viewpoint  above  leading  to  (89)  is  precisely  that  in  th in- 
wire theory  from  which  the  usual  approximate  kernel  stems,  and, 
indeed,  for  a constant-radius  wire,  (39)  reduces  to  the  usual 
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kernel.  Selected  sample  calculations  employing  (89)  yield 
results  which  exhibit  no  significant  error  compared  with  re- 
sults based  upon  the  exact  kernel  (83),  and,  therefore,  in  the 
interest  of  e f f i c iency, ca leu la t ions  reported  here  are  based 
upon  ( 89 ) . 

f.  S t ep p ed-Rad ius  Antenna 

If  one  wishes  to  determine  the  current  on  a stepped- 

radius  antenna,  driven  at  z=z  by  an  emf  source  of  V volts. 

g 

he  simply  replaces  of  (88)  by  V6(z-z  ) as  suggested  in 

z g 

(69).  With  this  replacement,  the  linear  system  (87)  can  be 
solved  for  the  antenna  current. 

g.  Scatterer  or  Antenna  above  Ground 

To  analyze  a s t epped-r ad ius  structure  above  and  perpendic- 
ular to  a perfectly  conducting  ground  screen  as  depicted  in 
Fig.  12,  one  only  has  to  appeal  to  symmetry  and  to  modify  the 
kernel,  (83)  or  (89),  in  the  usual  way. 

3.  SAMPLE  RESULTS 

In  Figs.  13-24  are  presented  selected  results  of  currents 
on  stepped-radius  wires.  In  all  cases  the  presence  of  the 
step  influences  the  magnitude  of  current  relative  to  wt;at  it 
would  have  been  on  a constant-radius  wire,  but  the  effect  on 
the  shape  of  the  distribution  depends  upon  length  and  certain 
features  of  the  excitation. 

Near  resonance,  where  the  forced  response  is  dominated 
by  the  resonant  response,  the  distribution  is  little  affected 
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by  Che  presence  of  the  step  or  its  position  on  Che  wire.  This 
is  supported  by  the  data  of  Figs.  13-15  where  are  displayed 
currents  on  wires  which  have  steps  in  various  locations  and 
which  are  illuminated  by  a normally  incident  plane  wave. 

Those  knowledgeable  in  thin-wire  theory  are  aware  that  the 
forced  current  is  dominant  on  a IX  wire  which  is  symmetrically 
excited  and  they  know  that  this  response  is  essentially  a 
shifted  cosine  function.  Only,  if  the  forcing  function 
possesses  an  odd-function  component,  can  the  resonant  cur- 
rent respond  on  a IX  wire  of  constant  radius.  On  the  other 
hand,  in  the  case  of  a IX  s tepped-radius  wire  having,  of 
course,  the  odd  symmetry  within  the  structure  itself,  one 
does  expect  a component  of  resonant  current.  This  is  vividly 
exhibited  by  the  data  of  Fig.  16  where  one  sees  a significant 
change  in  the  real  current  on  a s t epped-r ad ius  wire  over  that 
on  a constant-radius  wire. 

Fig.  17  shows  the  current  on  a s tepped-radius  scatterer 
caused  by  plane  wave  illumination  impinging  upon  the  wire  at 
angles  0 (w.r.t.  z axis).  Figs.  18  and  19  provide  infor- 

mation on  the  variation  of  peak  current  on  s t epped- r ad ius 
scatterers  of  lengths  below  and  above  X/2. 

Figs.  20-23  give  comparative  data  for  IX  antennas  excited 
symmetrically  and  antisymmetrically.  Again  one  observes 
little  influence  on  current  due  to  the  presence  of  the  step 
at  s=0  when  the  excitation  is  an  odd  function, but  a marked 
effect  is  seen  in  the  case  of  even- f unc t i on  excitation. 
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Lastly,  shown  in  Fig.  24  is  the  current  on  a stepped- 


radius  monopole  above  a ground  plane.  The 
of  such  an  antenna  suggests  that  no  major 
of  current  is  to  be  expected,  and  this  is 
results  of  Fig.  24. 


inherent  symmetry 
change  in  shape 
corroborated  by  the 
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iSiire  13.  Current  Induced  on  a Half-Wavelength,  S tepped-Rad  ius 
Scattercr  by  a Normally  Incident  Plane  Wave 


Figure  14.  Current  Induced  on  a 4X/10  S t epped-Rad i us  Scatterer 
by  a Normally  Incident  Plane  Wave 


r 


Figure  15.  Current  Induced  on  7X/10  Stepped-Radius  Scatterer  by  a 
Normally  Incident  Plane  Nave 


Figure  16.  Current  induced  on  a One-Wavelength  Stepped-Radius 

and  a Constant-Radius  Scatterer  by  a Normally  Incident 
Plane  Wave 


Figure  20.  Curretic  on  Symmetrically  Excited,  One -Wavelength  Uniform- 
Radius  A n t e n n a 


Kiglire  21.  Current  on  Symme t r icn  1 1 y Excited,  One-Wavelength,  Stepped 
Rad i nu  Antenna 


i’igure  22.  Current  t)n  A n t i - S ymrae  t r i <■  a 1 1 y Excited,  One-Wavelength 
Uniform  Radius  Antenna 
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SECTION  IV 


AN  ITERATIVE-SOLUTION  ANALYSIS  OF  THE  BENT-WIRE  SCATTERER 

1.  INTRODUCTION 

In  this  Section,  a theory  is  developed  by  means  by 
which  one  may  determine  the  current  induced  on  a bent-wire 
scatterer  illuminated  by  a specified  incident  electromagnetic 
field  (E^,H^).  The  solution  procedure  is  of  the  iterative 
type  and,  in  general,  follows  the  principles  set  forth  in 
recent  work  by  King  and  Wu  [8,14].  However  there  is  a sig- 
nificant difference  in  the  iterative  procedure  applied  to 
the  bent- wire  problem  compared  with  that  applied  to  either 
a straight  wire  or  perpendicular  crossed  wires.  This  analysis 
is  applied  to  the  problem  of  a bent-wire  scatterer  in  free 
space,  as  depicted  in  Fig.  25,  as  well  as  that  of  a bent- 
wire  scatterer  above  an  ideal  ground  plane  (Fig.  26). 

Calculated  results  are  presented  for  several  cases  of  interest. 

The  straight  wire  elements  of  the  bent-wire  structure 
are  treated  as  perfect  conductors  and  the  usual  simplifi- 
cations of  thin-wire  theory  are  employed  in  the  analysis. 
Principally,  the  element  radii  are  assumed  to  be  very  small 
compared  with  the  wavelength  a,  and  the  current  density  is 
considered  to  be  essentially  uniform  around  the  periphery  of 
an  element . 
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In  the  analytical  solution  procedure  which  follows,  an 

expression,  suitable  for  iteration,  is  derived  for  the  total 

axial  current  on  each  element.  In  Fig.  25  , one  sees  the  bent 

wire  located  in  a Cartesian  coordinate  system  with  the  bend 

at  the  origin.  One  wire  is  coaxial  with  the  z axis,  and  the 

other  extends  from  the  origin  radially  outward  in  the  yz 

plane  along  a line  which  is  declined  from  the  z axis  by  an 

angle  6.  Tne  direction  of  the  axis  of  the  latter  wire  is 

designated  the  s direction.  I^  is  the  total  axial  current 

on  the  z-directed  wire,  and  I is  that  on  the  s-directed 

s 

wire  as  shown  in  Fig.  25  which  serves  to  define  other  geom- 
etric quantities  of  interest. 


2.  ANALYSIS 


In  the  present  analysis,  expressions  for  I^  and  I^ 

are  formulated  and  are  used  to  obtain  zero  and  first  order 
solutions  by  iteration.  The  above-mentioned  expression  for 
I^  is  developed  below  from  consideration  of  the  magnetic 
vector  potential  A and  the  partial  differential  equation  which 
this  quantity  satisfies. 

Since  the  current  on  the  bent  wire  is  in  the  yz  plane, 
one  can  readily  show  [5]  that 


+ k‘ 


5_ 
3 z 


(90  ) 


where  A and  A are  the  z and 

z y 


y components  of  the  vector 


potential,  is  the  z component  of  the  electric  field  result- 

ing from  currents  induced  on  the  scatterer,  oj  is  the  angular 
frequency  of  the  suppressed  harmonic  variation  in  time,  and 

k is  2tt/X.  Consider  now  that  E ,A  , and  A are  evaluated  on 

z z y 

the  surfaces  of  the  z-directed  wire  in  Fig.  25  along  the 
line  (a^,0,z)  for  zs(-h,0).  In  order  to  satisfy  the  boundary 
condition  that  total  electric  field  tangential  to,  and  eval- 
uated on,  the  perfectly  conducting  wire  surface,  one  requires 
Ez ( 3z , 0 , z ) +E^ ( az , 0 , z ) =0  , ze(-h,0)  where  E^  is  the  z component 

of  E^.  Equation  (90)  is  a well-known  inhomogeneous  differ- 
ential equation,  and  one  can  demonstrate  [ 5 ] that,  on  the 
7-directed  wire, 

A (a  ,0,z)  = C cos  kz  + B sin  kz  + V^(z)  + U (z)  (91) 

Z Z ^ 7T  Z Z Z Z 


90 


The  2-directed  vector  potential  in  (91)  is 


A (a  , 0 , z ) = 7— 

2 Z 4TT 


/ 


I^(z ' }K(z-z ' ,a^)dz  ' 


z ' =-h 


+ COS0  

4tt 


/■ 


I^(s')G(z,s' ,a^)ds' 


s '=0 


where 


K(C,a) 


-jk[C^+a^ 


. 2 
[?‘+a  . 


'2 


and 


G(5,C  ,a) 


^-jk[a^  + 5'^sin^6  +(?-C'cos9)'^] 


[ a + C'~sin^9  +(5-l'cos0)^] 


At  a point  in  space  (a^,y,z). 


A is 

y 


Ay ( X , y , 2 ) 


I 


'4 

4tt 


s in0 


/ 


= 0 


(93) 


(94a) 


(94b) 


(95) 


91 


where 


R = [a“  + (y-s'  sinQ)“  + (z-s'  cos0)^] 

y ^ 


(96 


The  derivative  of  appearing  in  (92b)  is  seen  to  be 


37 


y = 0 


sin“0  f s ’ 
4tt  I 


I (s’), 
s 


s ' =0 


*3  ^ j 

R 


-jkR 


y ^ 


) ds  ' 

U = 0 


(97) 


Hence,  subject  to  interchange  in  order  of  integration  in  (97), 

U of  (92b)  becomes 
z 


2 

■sin  6 


s ' =0 


:=o 


?S 


k 

J 

Cs 


1 -jkR, 


cos  k(z-;)d;ds' 


(98) 


with 


i"  2 . ^ 

a^  + s sin“9  + (C“s'  cos3) 


^;s 


1, 

-il  i 


(99) 


?he  double  integral  of  (98)  can  be  reduced  [15]  to  the  single 


integral 
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i: 

r 

i 


I 

f. 

i 


i 


/■ 


- j kR 


Zs 


■sin"0  I s ' I (s') 

i f-, 

s'=0  ""  " 


® * R,  [a"  + s ' ^sln^9  ] 


( i;-s  ' cos  9 ) cos  k(z-?) 


•jR^^sin  k(z-?) 


ds  ' 


(100) 


i;  = o 


which  enables  one  to  express  in  the  font  of  (101),  which  is 


(z)  = -sin“9  I (s  ' ) I 2 


-/■ 


s '=0 


5 5 < ( z-s ' cos6 ) G (z  , s ' , a ) 

[a"  + s'“sin'^9]/|_  ^ 


0 9 

- j k [ a"  + s ' ■ ] 1 

+ s ' cos0K(s  ' ,a^)cos  kz  + j e " sin  kz  ) ds' 

(101) 

computationally  much  more  manageable  than  is  (92b).  At  z=0, 

U is  seen  from  (92b)  to  be  zero  and,  for  z less  than  zero  by 
z 

a few  radii  (a^^’  integrands  of  (101)  are  quite  easy  to 

integrate . 

With  the  kernel  K defined  in  (94a)  partitioned  into  real 
and  imaginary  parts, 

K(^,a)  = K^(f,a)  + jK.(;,a)  (102) 


where 
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K^(C.a) 


and 


K^(C.a) 


2 


sin  kv£  + a 


/T2~  2 

V 5 + a 


Equation  (91)  can  be  rearranged: 


/ 


z'=-h 


)K^(z- 


,a^)dz 


C cos  kz  + B sin  kz  + 
z z z 


U,„(z) 
z s 


U^j(3) 


where 


/■ 


U^g(z)  = -COS0  I Ig (s ' ) G(z ,s ' ,a^) ds ’ 


s'-O 


and 


U (z)  = -j  /*i  (z  ' )K.  (z-z  ' ,a  )dz  ' 


(103a) 


(103b) 


(z)  + 


(104) 


(105a) 


(105b) 


V — ' - ■■  ■ 

Now,  without  approximation,  (104)  yields 

f (z)I  (z)  = C cos  kz  + B sin  kz  + V^(z) 
z z z z z 

+ U^(z)  + U^g(z)  + U^^(z)  (106) 

subject  to  the  definition, 

0 

’i'z(z)  = 1-^  J*ij,(z')K^(z-z’ ,a^)dz*  (107) 

z ' =-h 

A procedure  parallel  to  that  leading  to  (104)  for  the  z- 
directed  wire  can  be  applied  to  the  s-directed  wire  to  obtain 

y (s)I  (s)  = C cos  ks  + B sin  ks  + V^(s)  + U (s) 
i s s s s s s 


+ u 

s z 


(s)  + (s) 


(108) 


where  and  B^  are  arbitrary  constants  peculiar 
directed  wire  and  where  Che  remaining  terms  in 
defined  below. 


to  the  s- 
(108)  are 
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s J 


(109a) 


V^(s) 


— [e^Cb  ,0, 

n J s'  s’ 


C)sin  k(s-C)dC 


(109b) 


U,^(s) 
s z 


•COS0  J l^(z* )G(s ,a’ ,ag)dz' 
z'=-h 


(109c) 


Usi(s)  = -j  J*ig(3 ’ )K^(s-s' ,ag)ds' 


(109d) 


s'-O 


U (s)  = -sin^e  nz(z')(  2 ^~~2  2 

® J V [a^  + z ' '^sin  6] 

z'=-h  ® 


(s-z'cos9)G(s,z*  »a^) 


-jk[a2  + z'2]  1 

+ z' co30K(z'  ,ag)cos  ks  + je  ® sin  ks  ) dz' 
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(109e) 


I 


a.  '^-Functions 

When  a<<X,  which  is  true  for  thin  wires,  the  real  part 
of  the  kernel  of  (94a)  is  highly  peaked  for  small  ^/\  and 
decays  very  rapidly  for  a.  These  properties  are  evident 
from  the  approximations  below: 

+ a^  ^ cos  + (a/X)^ 

+ xJi^lxA  + (a/X)^ 


Kj.(5.a)  = 


cos  k 


K^(C.a)  i 


J(Z/xA  + (a/X)2 


- I iT^y(C/X)2  + (a/X)^  . 


C/X  <<1  and  a/X  <<1 


(110) 


and 


^r(C.a)  = — , a/X  <<1  and  ?>>a 


(111) 


Due  to  the  sharply  peaked  behavior  of  near  5=0  and  the 
fast  decay  for  ^>>a,  exhibits  the  selecting  property  of  a 
delta  function.  The  ^-functions  defined  in  (107)  and  (109a) 
may  be  represented  by 
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The  'V- function  is  almost  constant  for  w remote,  in  terms  of 
numbers  of  a,  from  end  points  and  but  varies  quite  rap- 
idly within  a few  radii  of  either  end;  its  value  at  w^^  and  W2 
is  approximately  one-half  its  value  at  ^)  / ‘2.  for  (w2-Wj^)>>a. 

One  may  obtain  a good  approximation  of  'F'(w)  at  w^^  and  W2 : 


99 


a 


at 


and 


at  w=w^ 


For  reasonable  values  of  (w2~w^)>>a  but  not  exceeding  roughly 

2 

100,  4'(w)  is  not  significantly  greater  than  2 ^n  (~)  • Hence, 
for  f'(w)  typical  of  the  derivative  of  current  on  a thin  wire-- 
not  at  a free  end--the  second  term  in  (1^3)  ■<  s very  small  com- 
pared with  the  first  at  or  W2 , and  one  has  a simple  but 
good  approximation  of  t ' at  w^  and  W2: 

'T  (w, ) i - (114a) 

1 a 

r (wO  4 - i (114b) 

2'  a 

It  is  worth  noting  that  the  value  of  'i'  or  4"  at  either  w^  or 
W2  is  highly  sensitive  to  a but  is  insensitive  to  (w2~Wj^)  so 
long  as  (w2“Wj)>>a. 

In  view  of  the  properties  of  the  'i'-function  given  above, 

f (z)  and  4'  (s)  are  seen  to  be  almost  constant  except  near 
z s 

the  upper  and  lower  ends  of  the  individual  wires.  The  values 
of  4'^  and  4*^  may  be  obtained  from  (112)  with  appropriate 


interpretation  of  w ,w». 


and  a,  and  the  derivatives  of  these 


^-functions  at  the  wire  junction  (bend)  are,  from  (114), 


4'*  (0)  = - — (115a) 

z a 

z 

4'’(0)  - ^ (115b) 

s a 

s 

b.  Outline  of  Iterative  Procedure 

The  large  value  of  the  '{'-functions  for  thin  wires  and  the 

observation  that  some  terms  in  (106)  and  (108)  are  far  more 

significant  than  others  suggest  that  these  expressions  may  be 

amenable  to  an  iterative  solution  procedure  for  determining 

and  In  such  a procedure,  one  identifies  a so-called 

zero-order  solution  which  comprises  only  terms  that  themselves 

are  not  direct  functions  of  the  currents. 

The  zero-order  currents,  1°  and  1°,  from  (106)  and  (108) 

z s 

are 


’{'  ^z 

rt  ^ 


cos  kz  + 


B sin 
z 


kz  + 


V^(z) 


(116a) 


I°(s) 

s 


— \ C° 
'V  [ ^s 


o i 

cos  ks  + B sin  ks  + V (s) 
s s 


(116b) 


where  the  '{'-functions,  defined  in  (107)  and  (109a),  are 
approximated  by  (112).  One  observes  that  in  (106)  and  (108) 
the  U-functions,  each  of  which  depends  upon  current  on  the  bent 
wire  as  is  seen  explicitly  in  (101),  (105),  and  (109),  are 


! 

] 


101 


ignored  in  the  zero-order  solutions.  Expressions  (116)  are 

subjected  to  boundary  conditons  and  the  zero-order  coefficients, 

C°,  B°,  C°,  and  B°,  are  thereby  determined, 
z z s s 

Finally,  returning  to  (106)  and  (108)  and  using  the  zero- 

order  currents  as  approximations  to  the  exact  currents  in  the 

integrands  of  the  U-functions,  one  forces  and  to  satisfy 

boundary  conditions  and  thereby  calculates  C^,  B^,  , and 

Bg.  With  these  coefficients  in  (106)  and  (108),  one  has  what 

are  called  first-order  solutions  for  I and  I . 

z s 

In  a procedure  like  that  described  above,  the  influence 
of  the  zero-order  currents  upon  the  first-order  is  proportional 
to  the  factor  l/*f;  hence,  the  importance  of  a large  4'-function 
is  clearly  evident. 

The  boundary  conditions  which  one  enforces  are  (1)  that 
the  current  vanish  at  the  two  free  ends  of  the  bent  wire  and 
(2)  that  both  the  current  and  its  first  derivative  be  con- 
tinuous at  the  bend.  In  summary,  these  boundary  conditions 
are 


I,(-h)  = 0 


(117a) 


Is(^)  = 


0 


(117b) 


s = 0 


(117d) 


dz 


I,(z) 


z»0 


i;(o) 


Is(s) 


“ i;(o) 


c.  Current  and  its  Derivative  at  the  Bend 

If  the  structure  under  study  here  were  a straight  wire 
(bent  wire  with  0=0°)  or  even  the  complex  structure  compris- 
ing two  wires  which  cross  at  90°  [8, 14],  it  would  be  found  that 
an  iterative  solution  could  be  formulated  in  such  a way  that 
the  '{'-functions  would  be  essentially  constant  even  at  the 
confluence  of  wires.  However,  in  sharp  contrast  to  analyses 
of  the  structures  above,  the  iterative  solution  method  for 
the  bent  wire  must  fully  account  for  effects  of  the  large 
derivative  of  the  '{'-functions  at  the  bend.  In  addition,  and 
of  paramount  bearing  on  the  present  analysis,  one  cannot 
relegate  to  a secondary  role  the  contributions  from  one  wire 
to  the  value  of  current  on  the  other  when  the  point  of  eval- 
uation is  at  the  bend--a  point  common  to  both  wires.  In 

other  words,  the  determination  of  I and  its  derivative  at 

z 

z=0  cannot  be  made  without  cognizance  of  1^  at  s=0.  Hence, 
before  the  first  order  coefficients  are  calculated  from 
boundary  conditions  (117),  a thorough  investigation  of  cur- 
rent and  its  derivative  is  in  order. 

From  (106)  '{'  (0)1  (0)  is 
z z 

'{',(0)1, (0)  = C + U (0)  + U ,(0)  (118) 

z z z z s z i 

since  U^(0)  is  seen  from  (92b)  to  be  zero.  Notice  from  (105a), 


(94b),  and  (94a)  that 


is 


zs 


I 


U (0)  = -COS0  I I (s ' )K(s ' ,a  )ds ' 


s ' =0 


and,  even  though  the  integration  above  is  over  the  s-directed 
wire,  it  reduces  to 


U^g(O)  = -COS0  4'^(0)Ig(0) 


COS0 


(s' )K^(s 


s'=0 


a ) ds  ' 
z 


(119) 


Since  I (0)=I  (0)  as  required  by  (117c),  (119)  and  (118) 

z s 

yield 


'V  (0)[1  + cos0]I  (0)  = C -j  COS0  I I^(s')K^(s’  ,a^)ds' 
z z z 


/■ 


s'=0 


-j  J"  I^(z')K^(z'  ,a^)dz' 


(120) 


z ' =-h 


By  a similar  procedure,  one  obtains  the  following  for 

(0)1  (0)  : 

s ^ 


104 


0 


(0)[1  + cos9]I  (0)  = C -j  cos9  /I  (z')K  (z',a  )dz' 

S S S I Z X s 


/■ 


2 * *-h 


-j  (s ' )K^(s ' ,ag)ds ' 


s ' =0 


Subject  to  the  definitions  of  the  constants. 


P =~j  COS0 
z 


J*ig(s')K.(s',a^)ds'  + U^. 


(0) 


s ' =0 


and 


/■ 


•j  COS0  I ( z ' ) (z ' , a^ ) dz ' + U^^(O)  , 


z'=-h 


which  are  clearly  not  zero-order  contributors,  (120)  and 
lead  to 


C + P 

I (0)  = 5 

z''^-’  'i'  (0)  [1  + COS0] 


and 


C + P 

I (0)  = ^ 

(0)  [1  + cos9: 


(121) 


(122a) 


(122b) 


( 121  ) 


(123a) 


(123b) 


105 


1 

Next,  in  order  that  boundary  condition  (117d)  may  be 

enforced,  it  is  necessary  to  investigate  the  derivatives  of 

each  term  on  the  right-hand  sides  of  (106)  and  (108).  Only 

and  of  (106)  and  and  (108)  require  special 

attention:  nence  the  derivative  of  U (z)+U  (z)  at  z=0  is 

z zs 

investigated  in  detail. 

From  the  definition  of  U^(z)  in  (92b),  the  desired  deriv- 
ative of  this  function  is  seen  to  be 


= - s i n 6 — r K ( s ' , a ) 
os  z 


(0)  becomes 


106 


u;(0) 


/' 


sin  0 / Ig(s' ) K(s'  ,a^)ds’ 

s'=0 


(124) 


At  this  point,  attention  is  turned  to  the  derivative  of 

U (z),  which  can  be  written 
zs 


u;s(o)  = 


^ U (z) 
dz  zs 


] -co.e  [1^ 


G (z ,s ' , a ) 

-I  z = 0 


ds  ' 


s'=0 


One  can  show  that 


I7  G(z,s' ,a^) 


z = 0 


L 3s' 


= -COS0  -5 — - G(z,s',a  ) 


z = 0 


which  leads  to 


U'  (0) 

ZS 


/■ 


s'=0 


cos^0  I I^(s')  1^  G (0  , s ' , a^)  ds  ' 


Recognizing  that  G(0,s',a  )»K(s',a  ),  one  obtains 

z z 


u;s(o) 


COS  9 


s’ ) l—T  K(s' ,a  )ds' 

d S Z 


s ' =0 


(125) 


Now,  the  sum  of  (124)  and  (125)  is  simply 


U' (0)  + U'  (0) 
z z s 


K(s'  ,a^)ds' 


The  portion  of  the  integrand  of  (126)  involving  , the  real 
part  of  K,  can  be  integrated  by  parts  and  (126)  can  thereby 
be  converted  to 


cos  kH 


cos  ka 


o;(0)  o;^(o)  . i,(«)  1^(0)  — 


- 

J ds' 


I (s'  ) (s'  ,a,)ds' 

s r z 


s '=0 


■"J  / ^<=’>  It-  K^<s',a^)ds' 


in  which  the  excellent  approximation  K^(J,,a^)=  £°.L_JsA  is 

utilized.  Now,  in  view  of  the  selecting  property  of  the  f- 

cos  ka  ^ 

function  and  since r = — = -'i''(0)  as  indicated  in  (114), 

a ci  ^ 

Z Z 


one  arrives  at 


r 


u;(o)  + u;^(0)  i ig(^)  + ig(o)'^;(o)  - i;(o)4'^(o) 


/- 


s ' -0 


K^(s'  ,a^)ds' 


With  (127)  available,  the  derivative  of  (106)  is  readily 
to  be 


)l,=o  - [ 


4'^(z)i;(z) 


r (z)i^(z) 


z»0 


or 


f (0)i;(o)  + (0)1  (0) 

z z z z 


kB  + U' (0)  + U'  (0)  + U' . (0) 

z Z Z 9 Z 1 


= kB^  + i^a)  + ig(o)'i';(o)  - r (o)4'^(o) 

^ 0 
K^(S’  .a^)ds' 

s'=0 


■"1  / It-  "i''’’* 

2 ' =-h 


(127 

seen 


) dz  ' 
(128 


Since  I (0)»I  (0)  and  I'(0)=I'(0),  the  above  reduces  to 


(129) 


24*  (O)I'(O)  =>  kB  + I (Z)  + W 

z z ^ ' z s Z z 


where 


=■  j / 


2 * »-h 


K^(z-.a^)d2'  +j  Jl^(s’)  fjr  K^(s'a^)ds' 


s'=0 


(130) 


Similarly, 


where 


24'  (O)I'(O)  = kB  - I,(-h)  + W 

S S S Z Q s 


(131) 


/■ 


"s  “ J JP'  Ki(z'.a^)dz’  +j  Ip-  K.(s*,a^)ds' 


z ' =-h 


s '=0 


(132) 


From  (123)  the  requirement  that  current  be  continuous  at 

the  bend  (117c)  is  equivalent  to  the  following  relationship 

between  C and  C : 
s z 


C 

s 


4-^(0)  ["'z 


+ P 


P 

s 


(133) 


Similarly,  continuity  of  derivative  (117d)  requires 


(134) 


B 

s 


1 

L^z  k 


cos  kit 
I 


1 

k 


W 

z 


+ 7 I (-h)  ^ W 

k 2 ' ' h k s 

Now  the  need  to  consider  the  derivative  of  the  functions 

is  evident.  If,  for  example,  'l'^(O)  were  ignored,  and,  hence, 

did  not  exist  on  the  left-hand  side  of  (128),  the  term 

ls(0)'l'^(0)  on  the  right-hand  side  of  (128)  would  have  remained 

and  would  appear  in  (134).  Similarly,  a term  involving 

I (0)4'’  (0)  would  appear  in  (134),  and  the  two  together  would 
z s 

contribute  a term  proportional  to 

1^(0)  [f;<0)  - . -1^(0)  • 

j z s ■' 


Clearly,  a very  large  term  like  the  above  would  incorrectly 
"uncouple"  the  coefficients  and  B^  and  would  cause  B^  to  be 
prohibitively  large  in  magnitude, 
d.  Zero-Order  Currents 

The  zero-order  coefficients  are  determined  from  the  ap- 
plication of  the  boundary  conditions  (117)  to  Eqs.  (116) 

which,  since  V^,  V^,  4—  V^,  and  4—  are  zero  at  Che  bend, 
z s dz  z ds  s 

results  in 


c°  cos  kh  - B°  sinkh  = -V^(-h)  (135a) 

z z z 


cos  k£  + sin  ki.  ■ -V  (JZ;) 
s s s 


^ c°  = ^ c° 

'V  f ‘'s 

z s 


4<  z " H*  ®s 
z s 


where  the  f-f unctions  are  treated  as  constants  calculated  from 
the  approximations  (112)  evaluated  at  the  midpoint  of  the  ap- 
propriate wire.  Equations  (135)  can  be  solved  for  the  zero- 


order  coefficients 


4' 

^ V^(X,)sin  kh  + V^(-h)sin  kil 

T S Z 

s 

sin  k(h+X.) 


^ V^(JI.)cos  kh  - V^(-h)cos  kJ. 

T S Z 

s 


sin  k(h+X.) 


where  (h+il)ynX/2,  n-1,2,...,  . When  C°  and  B°  are  determined, 

z z 

C°  and  B°  follow  from  (135c)  and  (135d)  or 
s s 


C°  - C° 

S T Z 

z 


s 4'  z 

z 
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'!  Equation  (135d)  follows  from  boundary  condition  (117d)  which 

, equates  the  derivatives  of  the  currents  at  the  bend.  In  the 

computation  of  these  derivatives  for  determination  of  the 
zero-order  coefficients,  the  '{'-functions  are  treated  as 
constants  even  though  this  is  clearly  not  the  case  near  the 
bend.  However,  Eq . (135)  is  exact  within  thin-wlre  approx- 

imations and  supports  the  treatment  of  as  a constant  even 
at  the  bend. 

e.  First-Order  Currents 

Returning  to  (106)  and  (108)  and  using  1°  and  1°  of  the 

z s 

zero-order  solution  in  the  U-functions  where  I and  I 

z s 

appear,  respectively,  one  has  first-order  currents.  These 
first-order  solutions  are  improvements  based  upon  the  zero- 
order  approximations.  The  first-order  coefficients  C^,  , 

C , and  are  determined  from  requiring  (106)  and  (108)  to  i 

satisfy  the  boundary  conditions  (117).  To  fully  account  for  j 

the  fact  that  the  '{'-functions  possess  appreciable  derivatives 

at  the  bend,  the  enforcement  of  (117c)  and  (117d)  is  achieved 

through  (133)  and  (134).  Equating  I (H)  and  I (-h)  of  (106)  j 

S Z .j 

and  (108)  to  zero  and  eliminating  and  B^  by  means  of  (133)  i 

and  (134),  one  finds  that  C and  B satisfy  i; 

z z 

.i 

1 

I 5 ) 

1 cos  kh  - B^  sin  kh  =>  -F^(-h)  (137a) 


r 


where 


r 1 

— (0~)  Pg  k ~ ^3^^^ 


-P,  COS  kJl  - — sin  k£ 
2 k 


F^(z)  = V^(z)  + U^(z)  + U^g(z)  + U^^(z) 


F^Cs)  - V‘(s)  + 0,(s)  + V^^(s)  + U^j(s) 


The  solutions  of  (137)  are 


-F  (-h)  sin  ki!-  + T sin  kh 
z 

sin  k(h+X.) 


T cos  kh  + F^(-h)  cos  kJ, 
sin  k(h+il) 


Knowing  and  from  (139)  and,  sub s equen t ly , and  from 
(133)  and  (134),  one  can  calculate  and  from  (106)  and 


(108)  . 
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3.  ZERO-ORDER  CURRENTS  FOR  BENT-WIRE 
SCATTERER  ABOVE  GROUND  PLANE 

The  currents  and  on  the  bent-wire  scatterer  above 
a ground  plane,  depicted  in  Fig.  26,  can  be  calculated  in  a 
way  similar  to  that  described  in  IV-(ld).  The  single  dif- 
ference is  that  the  boundary  condition  (117a)  is  replaced  by 

I^(-h)  =■  0 (140) 

Subject  to  (140)  and  ( 11 7b ) - ( 11 7d ) , the  coefficients  C°  and 

z 

B°  for  the  zero-order  currents  on  the  structure  of  Fig.  26 
z 

are  found  to  be 


I 

^ V^(£)  cos  kh  + (-h)  sin  kJl 

T s k z 

s 

cos  k(h+<l) 


(141a) 


f (-h)  cos  kil  + 
k z 


Vg(Jl)  sin  kh 


cos  k(h+Jl) 


(141b) 


The  coefficients  C°  and  B°  follow  from  (141)  in  (136)  and, 
with  the  four  coefficients  known,  one  may  calculate  the  zero- 
order  currents  for  the  bent  wire  above  a ground  plane  from 
(116). 
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RESULTS 


a.  Bent  Wire  in  Free  Space 

Shown  in  Fig.  27  are  curves  of  zero-order  current  on  the 
bent  wire  of  Fig.  25  for  various  values  of  bend  angle  9.  The 
radii  are  equal  as  are  the  lengths  of  the  two  wires 
(a^  = ag*0.  OOlX  , 2»h=*A/8).  The  plane-wave  incident  field  is 
z-directed,  the  direction  of  propagation  is  parallel  to  the 
X axis,  and  the  magnitude  of  the  incident  electric  field  is 
adjusted  so  that  eH  = 1 volt.  In  Fig.  28  is  displayed  current 
calculated  for  the  same  structure  and  illumination  by  means 
of  an  accurate  numerical  technique.  Attention  is  called  to 
the  good  agreement  even  though  the  data  of  Fig.  27  are  only 
zero-order  solutions. 

b.  Bent-Wire  above  a Ground  Plane 

In  Figs.  29  and  30  are  shown  zero-order  current  distrib- 
utions on  the  bent-wire  scatterer  of  Fig.  26.  The  radii  and 
lengths  are  equal  (a^^a^'O . 00 635X  and  H=h=X/2),  and  the  incident 
field  is  that  described  above.  Numerical  technique  solutions 
for  the  current  on  this  scatterer  are  provided  in  Fig.  31, 
where  again  one  observes  very  good  agreement.  The  flatness  of 
the  current  for  9-60°  is  remarkable  as  is  its  sensitivity  to 
small  changes  in  9 about  60°  (Figs.  30  and  31).  Notice  also 
in  Figs.  30  and  31,  for  the  9=90°  case,  that  the  top-loading 
element  (s-directed  wire)  bears  very  little  current--zero 
zero-order  current  (Fig.  29)--and,  hence,  is  an  ineffective 
top  load . 


116 


The  zero-order  solutions  presented  here  are  good  ap- 
proximations to  actual  currents  and  such  would  be  true  for 
numerous  cases.  However,  the  goodness  of  the  zero-order 
solution  is  highly  dependent  upon  wire  lengths.  One  should 
not  expect  the  zero-order  results  for  the  bent  wire  in  free 
space  to  be  acceptable  for  Jl+h  fairly  near  odd  multiples  of 
X/2  or  near  odd  multiples  of  X/4  in  the  case  of  the  bent 
wire  above  ground.  Of  course,  the  first-order  solutions  should 
provide  accurate  results  except  at  odd  multiples  of  X/2  in 
the  free  space  case  and  at  odd  multiples  of  X/4  for  the  struc- 
ture above  ground. 
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I'igure  27.  Zero-Ordor  Current  on  Bent-Wire  Scatterer  Illuminated 
by  Normally  Incidinit  I'lane  Wave 


Figure  28.  Current  on  Bent-Wire  Scatterer  Illuminated  by 
Normally  Incident  Plane  Wave 
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(31V3S  3AI1V13U)  lN3yynD 
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Figure  29.  Zero-Order  Current  on  Bent-Wire  Scatterer  above 

Ground  Plane  Illuminated  by  Normally  Incident  Plane 
Wave 


(31V3S  3AI1V13H)  ITOnO 
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Figure  30.  Zero-Order  Current  on  Bent-Wire  Above  Ground  for 
Small  Bend  Angles  about  6 = 60® 


Figure  31.  Current  on  Bent-Wire  Scatterer  above  Ground  Plane 
Illuminated  by  Normally  Incident  Plane  Wave 


SECTION  V 


MEASUREMENT  FACILITIES  AND  APPARATUS 

The  experimental  measurements  reported  here  consist  of 
current  and  charge  distributions  on  thin  wire  structures.  The 
first  set  of  measurements  were  made  on  scattering  elements, 
the  second  on  a driven  antenna.  Both  types  of  structure  utilize 
an  image  plane  geometry,  since  this  geometry  allows  complete 
separation  of  the  test  structure  and  the  instrumentation.  The 
facilities  and  apparatus  used  for  the  measurements  are  described 
in  this  Section. 

The  frequency  at  which  the  measurements  were  made  was 
determined  by  three  factors:  (1)  the  size  of  the  image  plane, 

(2)  the  requirement  that  the  test  structures  be  compatible 
with  thin  wire  theory,  and  (3)  the  physical  size  of  probes 
and  associated  apparatus  which  could  be  accurately  constructed. 

Requirements  2 and  3 above  indicated  that  a frequency 
range  from  300MHz  to  600MHz  was  desirable.  However,  the 
then  existing  12  ft.  x 13  ft.  image  plane  was  too  small  for 
measurements  in  this  frequency  range  and  an  enlargement  of 
this  facility  was  necessary. 

An  18  ft.  X 28  ft.  image  plane  was  constructed  by 
extending  the  original  plane  to  the  fullest  extent  possible 
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in  the  antenna  laboratory  building.  The  new  portion  was 
constructed  of  4 ft.  x 12  ft.  x 0.032  in.  aluminum  sheets 


! 

i 

t 


cemented  to  5/8  in.  thick  particle  board  backing.  The  joints 
between  the  aluminum  sheets  are  backed  by  aluminum  strips  and 
10/32  flat  head  aluminum  bolts  are  counter  sunk  along  the  edges 
at  a spacing  of  about  6 in.  The  joints  and  bolt  heads  are 
covered  with  0.002  in.  thick  silver  backed  tape.  The  final 
structure  deviates  from  a plane  by  less  than  0.02X  overall 
at  600MHz  and  less  than  0.007X  over  short  distances  at  the 
same  frequency.  Mounting  for  the  test  antennas  is  provided 
at  a point  8.2  ft.  from  the  lower  edge  and  10  ft.  from  one 
side.  A cross  sectional  view  of  the  antenna  laboratory  showing 
the  position  of  the  image  plane  is  provided  in  Fig.  32. 

Panels  of  24  in.  high  pyramidal  absorber  can  be  placed  around 
the  image  plane  to  reduce  reflections, 
b.  Instrumentation 

The  technique  used  for  measurement  of  current  and  charge 
distributions  was  basically  the  same  as  that  used  by  many 
researchers  (see  for  example  Mack  [16]).  Charge  density  was 
measured  by  a very  short  monopole  perpendicular  to  the  surface 
of  the  test  structure.  The  accuracy  to  which  the  response  of 
this  probe  corresponds  to  charge  density  depends  primarily  on 
the  constancy  of  the  normal  component  of  electric  field  over 
the  length  of  the  probe.  Errors  may  therefore  be  expected 
very  near  changes  in  the  surface  of  the  test  structure  such 
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FRONT  VIEW 


Figure  32.  Antenna  Laboratory  and  Image  Plane 
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as  bends,  junctions,  or  ends  of  the  tubular  structures.  The 
current  distributions  were  measured  by  a small,  center-loaded, 
image  plane  loop.  The  probe  is  mounted  with  its  plane  in 
an  axial  plane  of  the  tubular  antenna.  The  loop  thus  responds 
to  the  component  of  magnetic  field.  The  accuracy  to  which 
the  response  of  this  probe  corresponds  to  current  density 
depends  on  several  factors.  Since  the  electric  field  is 
essentially  perpendicular  to  the  surface,  the  amount  of  elec- 
tric field,  or  dipole  mode,  response  which  is  added  to  the 
magnetic  field  response,  is  dependent  on  the  accuracy  to  which 
the  load  is  centered  at  the  top  of  the  loop.  This  location 
of  the  load,  however,  cannot  discriminate  against  response 
to  tangential  components  of  electric  field  and,  hence,  errors 
similar  to  those  of  the  charge  probe  may  be  expected  near 
changes  in  the  surface  of  the  test  structure.  Finally,  the 
finite  size  of  the  loop  results  in  errors  at  points  where  the 
current  is  changing  rapidly  with  distance  along  the  surface. 


or  where  the  tangential  magnetic  field  changes  rapidly  in  a 
direction  perpendicular  to  the  surface.  Probes  of  the  type 


used  here  have  been  analyzed  by  Whiteside  [17].  The  design 
of  the  probes  and  probe  carriages  are  shown  in  Fig.  33. 

The  test  structures,  both  scatterers  and  driven  antennas, 
were  constructed  of  cylindrical  brass  tubing.  A slot  1/16  in. 
wide  was  cut  axially  in  this  tubing  to  allow  movement  of  the 
probes  along  the  structure.  The  probe  carriage  was  designed 
to  slide  inside  of  the  tubes  and  extend  through  the  slot  to 


the  Surface  of  the  tube.  The  charge  probe  was  constructed 
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on  Che  probe  carriage  from  0.034"  O.D.  semi-rigid  coaxial 
cable.  The  outer  conductor  of  the  cable  is  soldered  to  the 
probe  carriage  and  cut  at  the  surface.  The  inner  conductor 
continues  beyond  the  surface  to  form  the  monopole  antenna. 

The  current  probe  was  constructed  from  0.02  in.  semi-rigid 
coaxial  cable.  The  cable  was  first  formed  into  a semi  circle, 
and  then  the  center  conductor  was  removed  for  a distance  of 
about  0.020  in.  starting  from  a point  exactly  at  the  top  of 
the  loop.  The  cable  on  one  side  of  the  loop  is  short  cir- 
cuited, the  other  side  is  joined  to  the  signal  cable.  The 
load  for  the  loop  antenna  is  thus  the  input  impedance  of  the 
coaxial  cable. 

Signals  from  the  probes  are  carried  by  coaxial  cables 
passing  inside  of  the  tubular  test  structures  to  the  measuring 
Instruments  located  behind  the  image  plane.  A block  diagram 
of  the  system  is  shown  in  Fig.  34.  Preliminary  to  actual 
measurements,  the  entire  system  including  probes  was  tested 
in  a coaxial  line  having  a slotted  tubular  inner  conductor. 

The  commerical  instruments  used  in  the  system  are  listed  be- 
low : 

Power  Source  - General  Radio  Model  1209  B Oscillator 

Hughes  Model  1401  H Amplifier 

Boonton  Model  230  A Amplifier 

Measurement  - Hewlett  Packard  Model  8405A  Vector  Voltmeter 
Hewlett  Packard  Model  5345A  Counter 

Recording,  Processing  and  Plotting  of  Data-  Hewlett 

Packard  Model  9820  Calculator  with  Model 
9862A  Plotter 
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Figure  3A.  Block  Diagram  of  Current  and  Charge  Measurement  System 


The  primary  purpose  of  this  measurement  program  is  to 
obtain  data  for  comparison  with  theoretical  solutions  of 
wire  models  (Idealized)  of  aircraft.  In  the  cases  of  the 
scattering  measurements,  it  is  therefore  very  desirable  for 
the  illumination  to  correspond  as  near  as  possible  to  that 
of  a plane  wave.  For  the  structures  studied  here,  this  re- 
quirement reduces  to  having  the  component  of  electric  field, 
which  lies  in  the  plane  of  the  scatterer,  perpendicular  to 
the  image  plane  and  constant  in  both  magnitude  and  phase  over 
the  area  of  the  scatterer.  In  the  actual  experimental  system 
the  transmitting  antenna  was  only  12  wavelengths  away  from  the 
test  structure  at  600MHz.  The  spherical  wave  front  pro- 
duced by  a short  dipole  at  this  distance  creates  a 15°  phase 
change  along  a 1 wavelength  long  scatterer.  Calculations 
show  that  a collinear  array  formed  by  mounting  a one-wavelength 
element  above  the  image  plane  and  driving  it  at  its  center 
should  produce  a phase  deviation  of  only  5°  over  a one  wave- 
length long  test  scatterer.  The  transmitting  antenna  actually 
used  consisted  of  this  collinear  type  element  mounted  in  a 
90°  corner  reflector.  In  addition  to  the  direct  wavefront  of 
the  transmitting  antenna,  additional  distortion  of  the  illum- 
inating field  may  be  caused  by  reflections  from  the  edges 
of  the  image  plane  and  surrounding  objects.  As  a test  of  the 
illuminating  field  two  monopole  scatterers  were  measured, 
one  of  length  h = A,  and  the  other  of  length  h = 1.5X.  The 


comparison  of  the  current  and  charge  distributions  measured 
with  theoretical  results  of  Burton  and  King  1X81  are  shown  in 
Fig,  35.  The  agreement  is  seen  to  be  reasonably  good  for 
the  h * A case,  but  significant  deviation  is  seen  for  h » 1.5A. 
This  may  be  taken  as  evidence  of  near  plane  wave  illumination 
to  a distance  of  one  wavelength  from  the  image  plane.  For 
scattering  elements  greater  than  one  wavelength  long,  errors 
due  to  non-plane  wave  illumination  may  be  expected. 

Two  types  of  scatterers  were  measured,  the  bent  wire  and 
V cross  shown  in  Fig.  36.  Each  of  the  structures  was  con- 
structed from  brass  tubing  0.25  in.  in  diameter  with  a wall 
thickness  of  1/32  in.  The  bends  and  junctions  were  formed  by 
soldering  the  tubes  together  with  a miter  joint.  Since  the 
current  and  charge  distributions  were  measured  at  all  points 
except  near  the  joints,  it  was  necessary  to  design  a system 
to  control  the  probe  position  past  the  joints  as  shown  in 
Fig . 37.  A small  Teflon  pulley  was  mounted  inside  of  the  tube 

at  the  end.  A pre-stretched  and  treated  nylon  dial  cord  was 
attached  to  each  end  of  the  probe,  the  portion  attached  to 
the  forward  end  of  the  probe  passed  around  the  pulley  and 
back  through  the  tube  to  a control  mechanism  behind  the  ground 
plane.  The  coaxial  cable  carrying  the  signal  from  the  probe 
was  a 0.023  in.  O.D.  semi-rigid  cable  capable  of  bending 
around  the  joint  inside  of  the  tube.  Initially,  a small  Teflon 
piece  was  inserted  at  the  bend  to  protect  the  coaxial  cable. 
Later  it  was  found  that  a small  build-up  of  solder  at  this 
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b.  Bent  Wire  - With  Sphere  at  Junction 


Figure  36.  Bent  Wire  and  V Cross  Scattarers 
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point  provided  satisfactory  protection  for  the  cable.  The 
probe  carriage  could  not  pass  through  the  junction;  theretore 
It  was  necessary  to  insert  the  probe  from  the  image  plane 
end  and  measure  up  to  the  junction.  The  probe  was  then  re- 
moved and  inserted  in  the  outer  end  of  the  antenna,  cables  and 
control  lines  being  inserted  through  the  slot,  after  which 
measurements  were  made  from  the  end  back  toward  the  junction. 
This  measurement  procedure  was  somewhat  time  consuming  and 
it  was  necessary  to  insure  that  the  system  remained  constant 
during  each  set  of  measurements.  Although  frequency  drifts 
were  observed,  they  were  less  than  4 parts  in  10^. 

Measurements  of  the  current  and  charge  distributions 
were  made  on  a thin  cylindrical  monopole  antenna  containing 
a step  in  radius.  The  ratio  of  radii  of  the  two  parts  of 
the  antenna  was  originally  intended  to  be  about  2 to  1. 

The  diameter  of  the  smaller  portion  was  restricted  by  probe 
size  to  be  at  least  0.25  in.,  thus  the  larger  portion  was 
required  to  be  at  least  0.5  in.  in  diameter.  Availability 
of  materials  governed  the  final  choice  of  0.625  in.  for  the 
larger  diameter  and  0.25  in.  for  the  smaller  diameter.  In 
order  for  the  larger  radius  to  be  compatible  with  thin- 
wire  theory,  it  was  necessary  to  make  measurements  at  a 
frequency  of  300MHz.  Scattering  measurements  were  not  con- 
sidered to  be  reliable  at  this  frequency  with  the  image  plane 
size  available.  A coaxially  driven  monopole  was  therefore 
used  for  these  measurements. 
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The  construction  details  of  the  coaxial  line  antenna 


system  are  shown  in  Fig.  38.  Current  and  charge  measurements 
were  desired  along  both  sides  of  the  step  in  radius.  For 
correlation  of  the  data,  it  is  essential  that  the  same  probe 
system  be  used  throughout.  The  larger-rad ius  portion  of  the 
antenna  was  therefore  placed  at  the  image  plane  so  that  the 
probe  could  be  withdrawn  from  this  portion  and  inserted  into 
the  smaller-rad ius  tube  from  behind  the  image  plane.  The 
0.625  in.  tube  forms  both  the  larger-radius  part  of  the  mono- 
pole and  the  inner  conductor  of  the  coaxial  line.  The  tube 
is  slotted  axially  over  its  entire  length.  Measurements  can, 
therefore,  be  made  inside  of  the  coaxial  line  for  determination 
of  antenna  apparent  impedance  and  normalization  factors 
needed  for  determination  of  absolute  values  of  current  and 
charge  on  the  antenna.  The  0.25  in.  diameter  tube  is  mounted 
coaxially  through  the  larger  tube  and  extends  beyond  it  a 
distance  of  25  cm.  to  form  the  smaller-radius  portion  of  the 
antenna.  Provisions  were  made  for  adding  unslotted  sections 
of  0.25  in  tube  to  extend  the  length  of  the  smaller-radius 
portion  of  the  antenna.  The  outer  conductor  of  the  coaxial 
line  has  an  inside  diameter  of  1.020  in.,  and  the  space  be- 
tween the  two  conductors  is  filled  with  a high  density 
styrofoam,  = 1.12.  The  characteristic  impedance  of  this 
line  is,  therefore,  27.8  ohms.  The  coaxial  line  has  an  over- 
all length  of  1 meter  (IX  at  300MHz)  and  is  short  circuited 
at  the  end  opposite  the  antenna.  Signal  is  applied  to  the 
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Figure  38.  Stepped-Radius  Coaxial-Llne/Antenna  System  Diagram 
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I line  by  a T connection  at  a distance  of  25  cm.  from  the 

I short  circuit.  The  system  before  and  after  assembly  is  shown 

I in  Fig.  39. 

I 

The  probes  used  for  these  antenna  measurements  are 
identical  to  those  described  in  the  previous  section.  The 
probe  carriage  differs  however,  since  only  straight  line 
movement  of  the  probe  is  necessary.  A 0.125  in.  diameter 
' brass  tube  is  attached  to  the  rear  of  the  probe  carriage  and 

I 

i extends  beyond  the  rear  of  the  coaxial  line.  A pointer 

, attached  to  this  tube  provides  for  measurement  of  probe 

I 

' location.  The  0.020  in  O.D.  semi-rigid  coaxial  line  from  the 

! 

• probe  passes  through  this  0.125  in.  diameter  tube. 

[ The  instrumentation  for  the  antenna  system  was  the  same 

f 

( 

! as  that  shown  for  the  scatterer  measurements. 

i 
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SECTION  VI 


EXPERIMENTAL  MEASUREMENTS — STEPPED-RADIUS 

The  apparatus  described  in  Section  V was  used  to  obtain 
current  and  charge  distributions  on  a thin  s t epped -rad ius 
monopole.  The  measured  data  is  presented  in  this  Section. 

For  interpretation  of  the  data  the  calibration  and  normali- 
zation factors  which  were  used  are  of  considerable  importance 
and  these  will  be  discussed  first. 

The  coaxial  line/antenna  system  provides  for  continuous 
movement  of  the  probes  into  the  coaxial  line  portion.  Measure- 
ments within  the  coaxial  line  make  possible  an  absolute  cali- 
bration of  the  probe  measurement  relative  to  the  voltage 
across  the  line  at  the  image  plane. 

A calibration  factor,  N , for  the  charge  probe  is  deter- 

q 

mined  to  satisfy  the  equation. 

N Q (z)  = q(z)  covlombs/meter-volt  (142) 

<\  n 

where  Q (z)  is  the  meter  reading  of  the  charge  probe  signal, 
m 

and  q(z)  is  the  linear  density  of  charge  on  the  antenna. 

Within  the  coaxial  line  the  charge  density  is  related  to  the 
voltage  by 

q(z)=*V(z)C  (143) 

o 
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where 
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Xn(b2  /b^) 
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(144) 


is  the  capacitance  per  length  of  the  coaxial  line.  The  most 
convenient  location  for  determining  q(z)  is  z = -A/2  where 


V(-A/2)  - 


From  Eqs.  (142)  and  (143)  we  then  have 


e ^0 


(145) 


q Q„(-^/2) 


(146) 


where  V(0)  = 1 volt  has  been  assumed.  The  capacitance  of  the 
coaxial  line  in  this  system  is  - 1.272  x 10~^®  farad/meter. 
The  calibration  factor  was  applied  to  all  measurements  on 
the  antenna. 

The  calibration  factor  for  the  current  probe  measure- 
ments was  obtained  in  a manner  similar  to  . At  any  point 
within  the  line 


I(z)  = V(z)Y(z)  (147) 

The  most  convenient  point  for  computing  N^.  is  z = -A/4,  where 
y (-1/ 4)  = 2(0),  and 


V(-X/4)Z(0)  , 

I (-A/4)(Z_,)2 
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(148) 
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where  I (-A/A)  is  the  meter  reading  of  the  current  probe 
m 

signal.  The  voltage  V(-A/4)  may  be  related  to  V(0)  by 


V(-X/4) 


2(0) 


V(0)e 


ji:/2 


(149) 


where  is  the  characteristic  impedance  of  the  coaxial 
line.  Combining  148  and  149 


^jTT/2 

“ I (-X/4)Z 
m c 


(150) 


where  V(0)  =•  1 volt  is  assumed.  The  calibration  factor 

determined  by  eq.  150  depends  only  on  knowledge  of  Z^. 

The  factor  could  be  computed  from  measurements  at  any  point 

within  the  line,  but  knowledge  of  the  antenna  impedance 

would  then  be  necessary.  The  disadvantage  of  eq.  150  is 

that  when  the  antenna  impedance  is  high,  as  occurs  for  one 

of  the  test  antennas,  the  measurement  of  I must  be  made 

m 

near  a current  minimum. 

The  probe  calibration  factors  discussed  above  provide 
for  absolute  normalization  of  measurements  on  that  portion 
of  the  antenna  which  is  a smooth  continuation  of  the  coaxial 
line.  Measurements  on  the  smaller  radius  portion  of  the 
antenna  require  additional  correction  factor.  If  it  were 
assumed  that  the  probes  respond  to  the  fields  at  the  surface 
of  the  antenna,  the  signal  from  the  probes  would  then  be 
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proportional  to  the  surface  density  of  current  or  charge. 


The  calibration  factors  N_  and  N convert  the  measurements 

I q 

to  total  current  and  linear  density  of  charge  respectively. 
Hence,  the  assumption  that  a probe  signal  is  proportional 
to  surface  density  leads  to  a correction  factor  which  is 
simply  the  ratio  of  the  two  radii,  a/b^^.  Unfortunately, 
the  situation  is  not  this  simple.  The  probes  being  of 
finite  dimensions  respond  to  fields  away  from  the  surface 
of  the  antenna,  and  this  must  be  considered  in  the  cor- 
rection factor . 

The  charge  probe,  being  a very  short  monopole,  responds 
to  the  integral  of  the  over  its  length.  At  points  along 
the  antenna  which  are  not  too  near  the  step  or  end,  the 
field  varies  1/r  near  the  surface  of  the  antenna.  Hence, 
the  charge  probe  response  is  proportional  to 
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dr 
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Jin[(R  + h)/R] 


(151) 


where  R = b,  or  a,  the  radius  of  the  antenna,  and  h = the 
length  of  the  charge  probe.  The  correction  factor,  M , to 

q 

be  applied  to  measurements  of  charge  on  the  small  radius 
portion  of  the  antenna  is  therefore 
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M 
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For  the  two  radii  used  » 0.469  with  a probe  height  h = 

0.1  in. 

The  correction  factor  for  the  current  probe  is  ob- 
tained in  a similar  manner.  The  probe  response  in  this 
case  is  proportional  to  the  integral  of  over  the  area 
of  the  loop  probe.  As  for  we  may  assume  that  varies 
as  1/r  near  the  surface  of  the  antenna  except  near  the 
step  and  end.  Thus  the  signal  from  the  probe  is  proportional 


ffi 


(153) 


For  evaluation  of  this  integral  the  coordinate  system  of 
Fig.  40  is  used  and  eq.  153  becomes 


’ 1 2 
h -(y-a)^ 


dxdy 


(154) 


Integrating  with  respect  to  x and  introducing  the  change 
of  variable,  y « z + R,  we  obtain 


F = 2 


1 2 , 
^ /h  + z dz 


(155) 


Expanding  the  factor  l/(z+R)  in  a povjer  series,  and  inte- 
grating term  by  term  a series  solution  is  obtained. 


Figure  40.  Coordinate  System  for  Loop  Probe  Corrections 


F . 2h  I (|)(1  t (A)")  - i(|)'  - (156) 

Taking  the  ratio  of  F(R»b)  to  R(R  = a)  and  us''ng  a probe 
radius  h * .0625  in.  the  correction  factor  is  found  to  be, 

Rj  = 0.446. 
c.  Measured  Data 

The  current  and  charge  distributions  for  two  stepped- 
radlus  monopoles  have  been  measured.  The  measurements  were 
conducted  at  a frequency  of  300MHz , and  the  electrical  radii 
at  this  frequency  are  kb^^  » 0.05  and  ka  = 0.02.  Fig.  41 
gives  the  data  for  h^^  = X/4,  h^  “ X/4.  Fig.  42  gives  the 
data  for  h^  = X/4,  h2  = X/2.  These  two  cases  were  choosen 
to  illustrate  the  behavior  for  charge  minimum  and  charge 
maximum  at  the  position  of  the  step.  Some  irregularities 
in  the  data,  particularly  the  charge  distribution  of  Fig.  42 
are  related  to  equipment  problems  which  had  not  been  over 
come  at  the  time  of  this  report.  The  behavior  of  the 
charge  data  near  the  step  shows  the  expected  edge  behavior 
of  exterior  and  interior  90°  corners. 

Probe  errors  are  significant  near  the  step  where  the 
fields  depart  from  the  simple  behavior  assumed  in  computing 
the  probe  response.  The  departure  of  the  current  distribu- 
tions from  a smooth  curve  near  the  step  is  believed  to  be 
entirely  caused  by  probe  error.  It  is  difficult  to  deter- 
mine at  tnta  time  the  degree  to  which  the  observed  charge 


beha'ijr  -^ur  the  step  is  influenced  by  probe  error.  The 
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measured  data  demonstrates  continuity  of  current  across  the 
step.  However,  until  the  probe  errors  can  be  determined 
no  definite  conclusion  can  be  reached  concerning  continuity 
of  charge  per  unit  length  at  the  step. 

Figs.  43  and  44  show  a comparison  between  the  measured 
current  and  the  numerical  solution  for  the  case  h^^  ■ h2  = X/4. 
The  agreement  is  seen  to  be  very  good  except  for  a small 
difference  which  may  be  the  result  of  errors  in  the  probe 
calibration . 
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Figure  43.  Stepped-Radius  Antenna,  hj^=X/4,  h2=X/4 
ll|.  Phase  I Theory  vs.  Experiment 
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Figure  44.  Stepped-Radius  Antenna,  h^=A/A,  n^=\/4 

I-,  I,.  Theory  vs.  Experiment 
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SECTION  VII 


EXPERIMENTAL  DATA-BENT  WIRE  AND  V-CROSS 

The  experimental  data  obtained  on  the  scattering  elements 
described  in  Section  V arepresented  in  this  Section.  The 
principal  use  of  thesedata  is  expected  to  be  comparison  with 
theoretical  solutions.  For  this  purpose  the  shape  of  the 
distribution  curves  is  the  most  important  Information.  No 
attempt  has  been  made  to  obtain  absolute  normalization  factors 
for  the  amplitude  or  phase  of  the  measured  data.  It  is, 
therefore,  suggested  that,  in  any  future  use  of  these  data, 
normalization  can  be  made  at  any  appropriate  point  on  the 
curves . 

The  geometries  and  coordinate  systems  used  for  presentation 
of  the  data  are  those  shown  in  Fig.  26  and  Fig.  45.  The 
actual  structures  upon  which  measurements  were  made  are 
shown  in  Fig.  36. 

The  measured  current  and  charge  distributions  for  bent 
wire  scatterers  are  given  in  Figs.  46,  47,  48,  49.  The 
current  distributions  for  the  case  h^  » h^  = ^/2  were  measured 
using  a frequency  of  600MHz,  at  which  the  electrical  radius, 
is  ka  * 0.04.  These  results  are  shown  in  Fig.  46,  and  it  is 
seen  that  the  current  on  the  9 = 30°  case  is  very  similar 
to  the  9=0,  i.e.  straight  wire,  case.  The  principal  difference 
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Figure  45.  Coordinate  System  for  V Cross  Scatterer 
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between  these  two  distributions  occurs  near  the  image  plane, 
together  with  a small  shift  in  the  position  of  the  maximum. 


Increase  of  the  bend  angle  to  0 = 62°  results  in  a very 

significant  change  in  the  distribution  of  current  on  the  z = 0 

to  -h^, portion.  It  is  somewhat  surprising  that  the  current 
distribution  on  the  z=»0  to  + h2  portion  retains  essentially 
the  same  shape  as  for  the  6 »»  0°  and  30°  cases.  A comparison 
with  Figs.  29  and  30  shows  very  good  agreement  with  the  pre- 
dictions of  the  zeroth  order  solution  as  well  as  the  numerically 
solutions  given  in  Fig.  31.  The  curves  shown  across  the  junction 
region  in  Fig.  46  are  of  course  only  estimates  of  the  behavior 
in  this  region.  In  addition  to  the  mechanical  problem  of  moving 
the  probes  into  the  junction  region,  measurement  errors 
become  significant  near  the  junction  where  the  fields  are 
not  slowly  varying  over  the  dimensions  of  the  probe.  These 

errors  are  evident  in  the  0 = 62°  data  of  Fig.  46  where  the 

two  data  points  nearest  the  junction  have  been  excluded  from 
the  estimated  curve. 

The  charge  distributions  given  in  Fig.  47  show  many  of 
the  same  features  as  discussed  above  for  the  current.  The  0*0° 
and  0 ■ 30°  distributions  are  seen  to  be  very  similar  with 
only  a shift  in  the  position  of  the  minimum.  The  0^62°  data, 
however,  shows  even  greater  differences  from  0 = 0°,  than 
were  evident  in  the  current  distribution.  In  fact  the  entire 
nature  of  the  charge  distribution  is  changed,  with  the  minimum 
which  occurs  at  the  junction  for  0=0°  and  30°  being 
replaced  by  a maximum., 
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Figs.  43  and  47  give  the  measured  current  and  charge 
distributions  for  bent  wires  with  These 

measurements  were  made  at  a frequency  of  300MHz  and  an 
electrical  radius  of  ka  » 0.02.  Although  not  shown  here  it 
may  be  seen  that  the  distributions  for  0 =30°  are  very  similar 
to  the  known  distributions  for  a straight  wire,  9-0°.  As  for 
the  longer  element  the  9=  62°  distributions  for  this  case 
show  significant  changes,  particularly  in  the  charge,  a 
contrast  between  the  h = X/2  and  h = X/4  charge  distributions 
for  0 = 62°  is  partiuclarly  evident  in  the  charge  minimum 
which  occurs  at  the  bend  for  h = A/4  whereas  the  h = A/2  case 
shows  a maximum  of  charge. 

All  of  the  bent  Wire  structures  discussed  above  were 
constructed  with  a sharp  mitered  junction.  A different 
junction  geometry  was  obtained  by  placing  a brass  sphere  of 
1 in.  diameter  around  the  junction.  It  must  be  pointed  out 
that  this  is  a rather  large  sphere  and  does  not  satisfy  the 
conditions  of  thin  wire  theory  at  the  frequencies  used.  The 
data  are  shown  in  Figs.  50  and  51  for  distributions  of  current 
and  charge  with  and  without  the  sphere.  Examination  of  these 
data  shows  that  for  the  case  of  charge  minor  changes  in  the 
distributions  occur.  However,  for  the  case  which  has  a charge 
maximum  at  the  junction  very  significant  changes  occur  when 
the  sphere  is  placed  over  the  junction.  For  junction  regions 
having  an  electrical  size  ka  = 0.16,  the  geometry  of  the  junction 
region  is  therefore  of  considerable  importance.  Further  in- 
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vestlgation  will  be  necessary  before  conclusions  can  be  reached 
concerning  junctions  of  small  electrical  size. 

The  V cross  structure  is  perhaps  the  most  important  in 
relation  to  theoretical  studies  of  aircraft.  The  experimental 

data  which  has  been  obtained  for  current  and  charge  distributions 
on  this  structure  is  given  in  Figs.  52,  53,  54,  and  55. 

These  measurements  were  made  over  an  image  plane  as  described 
in  Section  V.  The  two  cases  studied  wereintended  to  be  examples 
of  the  behavior  when  a charge  minimum  occured  at  the  junction, 
and  when  a charge  maximum  occured  at  the  junction.  As  seen 
in  Fig.  53  the  h^  = h^  “ X/2,  =X/4  case  achieved  a 

very  deep  charge  minimum  at  the  junction.  However,  as  seen 
in  Fig.  55  the  h^  =•  A/2,  h^^  “3A/4  case,  while  having  a significant 
charge  density  at  the  junction,  did  not  produce  a real  maximum 
desired.  The  two  cases,  however,  do  provide  examples  which 
show  significantly  different  behavior. 

Probe  response  errors  exist  near  the  junction,  but  do  not 
appear  to  be  significant  factors  in  the  measured  data.  A 
much  more  significant  error  is  believed  to  be  caused  by  the 
failure  to  achieve  a planar  wave  front  illumination.  It  is 
not  possible  at  this  time  to  evaluate  the  amount  of  error 
caused  by  the  unknown  incident  field;  however,  the  shape  of  the 
distribution  on  each  arm  is  believed  to  be  accurately 
represented  in  the  measured  data.  Some  changes  may  be  expected 
in  the  relative  magnitudes  of  these  currents  with  changes  in 
the  incident  field.  Additional  studies  of  these  structures 
j will  be  necessary  before  definite  conclusions  can  be  reached. 
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Figure  51.  Bent  '.fire  Scatterer  Current  and  Charge 

” 1*2  = X/2,  ka  = 0.04  with  and  without 

Sphere  at  Junction 
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Figure  54.  V Cross  Current,  = 3X/4,  h2  = >72 
“ ^2  “ ® “ 60° 
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SECTION  VIII 


CONCLUSIONS 


Results  of  the  theoretical  and  experimental  investiga- 
tions are  presented  for  current  and  charge  distributions  on 
thin-wire  scatterers.  Insofar  as  the  investigators  are 
aware,  these  are  the  first  experimental  data  available  for 
the  structures  studied  in  this  project.  The  study  of  wire 
junctions  is  far  from  complete  but  this  initial  work,  both 
theory  and  experiment,  certainly  addresses  the  major  question 
^ of  what  the  appropriate  thin-wire  junction  conditions  should 

be.  The  experimental  investigation  of  structures  with  junc- 
tions is  in  its  infancy  and  numerous  facets  deserve  further 
^ refinement,  yet,  within  limitations  of  accuracy  and  experi- 

mental error,  the  correlation  between  measured  and  theoretically 
[ predicted  results  is  both  gratifying  and  encouraging.  ; 

A careful  study  of  the  properties  of  currents  and  charges  , 

in  a coaxial  line  with  a step  in  its  inner  conductor  and  of  j 

i 

a spheroidal  antenna  with  a similar  step  leads  to  the  conclu-  I 

sion  that  the  behavior  of  the  charges  per  unit  length  near  i 

the  junction  of  two  cylindrical  conductors  with  different  j 


radii  is  not  simple  even  for  electrically  thin  conductors. 
However,  a good  approximation  is  given  by  the  condition  (See 


Eq.  (64),  p.  43) 
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Is  the  expansion  parameter  for  conductor  1 whose  definition 
is  somewhat  arbitrary,  but  always  involves  a logarithmic 
function. 

From  the  analysis  of  a stepped-radius  wire  scatterer 
presented  in  Section  III,  one  sees  that  the  presence  of  the 
step  influences  the  magnitude  of  current  relative  to  what 
it  would  have  been  on  a constant-radius  wire,  but  the  effect 
on  the  shape  of  the  distribution  depends  upon  length  and 
certain  features  of  the  excitation. 

Near  resonance,  where  the  forced  response  is  dominated 
by  the  resonant  response,  the  distribution  is  little  affected 
by  the  presence  of  the  step  or  its  position  on  the  wire. 

However,  on  a antiresonant  (IX)  wire  of  constant  radius  sub- 
ject to  even-function  excitation,  the  current  is  essentially 
the  shifted  cosine  forced  response.  Only,  if  the  forcing 
function  possesses  an  odd-function  component,  can  the  resonant 
current  respond  on  a IX  wire  of  constant  radius.  On  the  other 
hand,  in  the  case  of  a IX  stepped-radius  wire  having,  of 
course,  the  odd  symmetry  within  the  structure  itself,  a strong 
antiresonant  current  is  produced  by  an  e ven- f unc t ion  excitation. 

Measured  current  and  charge  on  a stepped-radius  struc- 
ture are  presented  in  Figs.  43  and  44  where  one  sees  good 
correlation  between  experiment  and  theory.  However,  measure- 
ment of  the  condition  on  charge  at  the  junction,  discussed 
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In  Section  II,  is  not  entirely  satisfactory.  The  investi- 
gators do  not  judge  the  accuracy  of  their  data  to  be  sufficient 
to  enable  them  to  claim  undeniable  corroboration  of  the  small 
discontinuity  predicted  in  linear  charge  density  for  all 
cases  considered. 

In  the  bent-wire  scatterer  above  a ground,  the  investi- 
gators found  far  greater  sensitivity  of  current  upon  bend 
angle  than  expected.  Actually,  in  an  attempt  to  uncover  rea- 
sons for  disagreement  between  experimentally  and  theoretically 
determined  currents,  a 2 degree  error  was  found  in  the  angle 
of  the  fabricated  bent  wire.  Correction  for  this  brought 
theory  and  experiment  into  good  agreement  (shape)  as  can  be 
seen  from  comparisons  of  Figs.  29-  31  with  Figs.  46  and  48. 

A small  conducting  sphere  (ka  = 0.16)  was  placed  at 
the  junction  of  the  bent-wire  structures  for  the  purpose  of 
observing  effects  of  junction  geometry.  The  measured  currents 
distribution  shows  little  change  when  a charge  minimum  occurs 
at  the  junction.  However  when  a charge  maximum  occurs  at 
the  junction  very  significant  changes  in  the  current  distri- 
bution are  observed. 
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